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A Complete Self-Contained Audio Measurement System, by James D. Foote This new 
audio analyzer has everything needed for audio measurements — source, filters, detectors, 
voltmeter, and counter — all under microprocessor control. 

Audio Analyzer Applications The major areas are general audio testing, trans- 
ceiver testing, and automatic systems. 

Making the Most of a Microprocessor-Based Instrument Controller, by Corydon J. Boyan 

fn an audio analyzer, microprocessor control means automatic operation, "guaranteed" 
accurate measurements, and extra features. 

Design for a Low-Distortion, Fast-Settling Source, by George D. Pontis it's based on a 

state-variable filter with refinements. 

Floating a Source Output, by George D. Pontis The fioating output lets the user eliminate 
ground-loop errors, sum signals, and add do offsets, 

A Digitally Tuned Notch Filter, by Chung Y. Lau ft eliminates the fundamental frequency 
component of the incoming signal for distortion and noise measurements. 

A Custom LSI Approach to a Personal Computer, by Todd R. Lynch Nine HP-produced 

large-scale integrated circuits make the HP-85 possible, 

Handheld Calculator Evaluates Integrals, by William M. Kahan Now you can carry in your 
pocket a powerful numerical integrator like those available on large computers. 

In this Issue: 

Audio frequencies are frequencies within the range of human hearing, roughly 20 Hz to 20 
kHz. Frequencies on either side of this range are often loosely classified as "audio, " too. 
Thus, the frequency range of the instrument shown on the cover, Model 8903A Audio 
Analyzer (page 3). is 20 Hz to 100 kHz, The 8903A is used for testing — among other 
things — many of the electronic devices that speak to us or play music, such as CB radios and 
high-fidelity systems. In our cover photo it's shown plotting the frequency response of the 
stereo amplifier on the left at different signal levels. 
There are, of course, other ways of making the measurements in the 8903As repertoire. 
What makes the 8903A better? First, its a complete system that includes a low-distortion signal source, a 
counter, a voltmeter, and various filters and detectors. Second, all of this is under microprocessor control, 
automatically stepping through complicated sequences of measurements and computations. Third, it's ex- 
tremely accurate; for example, it can measure total harmonic distortion (THD) down to 0.003% under normal 
conditions. Fourth, the 8903A has recorder outputs that make plotting results about as easy as it can be. 

ff you're not so fortunate as to have studied integral calculus in school, the integral of a function probably isn't 
the familiar and highly useful concept that it is to scientists and engineers. One way to think of an integral is this: 
draw the graph of the function as a meandering line on a piece of graph paper. Then the integral of the function is 
the area bounded by 1) the graph of the function, 2) the horizontal axis of the graph paper, and 3) two vertical 
lines called the upper and lower limits of integration. Sometimes the integral of the function can be expressed 
neatly in mathematical terms, but more often than not it can't. So various methods have been devised for 
estimating integrals using computers, Because most of these numerical integration programs run on very, very 
large computers, it seems like a small miracle that you can now carry a numerical integrator — a very good 
one — around in your pocket. Beginning on page 23, its designer tells us about its capabilities and limitations. 
On page 18 is an article about the nine special integrated circuits that make the HP-85 Computer possible. 
This set of custom integrated circuit chips minimizes the cost of the electronics, eases the problem of cooling the 
computer, makes the small package possible, and provides features that couldnt have been included otherwise. 

-R. P. Dolan 
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A Complete Self-Contained Audio 
Measurement System 

This automatic, autoranging audio analyzer has the signal 
source, distortion analyzer, and counter to make the 
measurements most often needed tn audio-frequency 
testing. 

by James D. Foote 



HEWLETT-PACKARD'S NEW MODEL 8903A Audio 
Analyzer (Fig, 1 1 is a complete audio measurement 
system for quick and accurate characterizati on 
of systems and signals in the frequency range 2CJ Hz to 100 
kHz, The starting point for the 8903 A is the classical distor- 
tion analyzer. Added to this are microprocessor control, a 
reciprocal frequency counter, rms detectors, and a pro- 
grammable audio SOUrce. These provide m < urate mea- 
surement of ac level, distortion, SIN AD/ signal-to-noise 
ratio, and dc level. The audio source and the measurement 
circuits can work independently or together. The source is 
programmable in frequency and level and has very low 
distortion. The measurement circuits can monitor this in- 
ternal source or any other independent input waveform. 
Together the source and measurement input can be used for 
swept response measurements. 
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All measurements are available at the push of a button 
No knob adjustment or operator interaction is necessary. 
One simply applies the signal and selects the measurement 
mode. All control and processing are handled by the inter- 
nal microprocessor. The microprocessor monitors the input 
signal and makes internal gain and frequency adjustments 
as required, 

In automatic measurement systems, the 8903 A is capable 
of rapid and straightforward remote control. Analyzer op- 
erations can be controlled and all measurements can be 
transferred via the Hewlett-Packard Interface Bus (HP-IB), 
Hewlett-Packard's implementation of IEEE Standard 4H8- 
1978. On the bench, the 8903 A allows rapid and accurate 
circuit characterization when many repetitive measure- 
ments are necessary. 

Major application areas for the BHOU A Audio Analyzer are 
general audio testing, transceiver testing, and automatic 
systems. In general audio testing, the 8903 A measures the 




Fig, 1. Model 8903 A Audio 
Analyzer makes the accurate 
measurements needed to charac- 
terize systems and signals m its 
frequency range of 20 Hz to 100 
kHz ft has applications in general 
audio testing, transceiver h 
and automatic systems Micro- 
processor control makes it auto- 
matic and easy to use. 
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Fig. 2. The 8903 A Audio Anafyzer is basically a distortion analyzer, with a tunable notch filter to 
remove the fundamental frequency component of the signal and a detector to measure what 
remains, which consists of noise and distortion. Added to this are a microprocessor-based 
controller, a reciprocal frequency counter, a programmable audio source, and other resources r 
forming a complete audio measurement system. 



frequency response and distortion characteristics of filters, 
high-quality amplifiers, audio integrated circuits, and other 
devices. The frequency of the internal source can be swept 
while making measurements in all modes, The analyzer 
provides recorder outputs and scaling for easy generation of 
plots using an X-Y recorder. 

For transceiver applications the most common receiver 
measurements are SINAD for FM receivers and signal-to- 
noise ratio for AM receivers. A psophometric filter is in- 
cluded for making measurements to CEPT standards. 
Common transmitter measurements such as audio distor- 
tion and squelch tones are made using the 8903A with its 
companion instrument, the 8901 A Modulation Analyzer, 1 
In automatic systems, the 8 903 A provides many frequently 
needed audio functions, doing the work of an audio syn- 
thesizer, digital multimeter, frequency counter, and tun- 
able notch filter. More details on specific applications are 
presented on page 6. 

Control Philosophy 

Front-panel control of the audio analyzer is simple, yet 
powerful. Most functions can be used and understood with 
very little training. The casual user can select amplitude. 
frequency, measurement mode, and filtering simply by 
reading the labels on the controls. More details are available 
on the instrument' s pull-out card. 

A great deal of measurement sophistication is built into 
the 8 903 A software. Measurement routines are structured 
to optimize measurement speed and accuracy. As a rule, 



measurements triggered from the bus or initiated from the 
keyboard are accurate from the first reading. The operator 
needn't wait for successive measurements to verify that the 
reading has stabilized. The software algorithm monitors 
key voltages in the audio chain and waits until they 
stabilize before taking data. Not only does the software 
perform these functions much more rapidly than the 
operator, but it can also ensure the optimal convergence of 
the measurement with a repcatable, well defined technique. 
Distortion, SINAD, and signal-to-noise ratio, in particular, 
are examples of measurements that in the past required a 
significant amount of settling time and operator interac- 
tion. A classical distortion analyzer requires repeated ad- 
justments to achieve an accurate distortion reading. More 
recent analyzers have offered semiautomatic tuning and 
leveling, but response time is often long, and operator in- 
teraction is required if the frequency, amplitude or relative 
distortion of the signal changes significantly. With the 
8903 A, some delay still exists, but the delay involved is 
minimized by careful circuit design and microprocessor 
control. 

Special functions extend user control of the instrument 
beyond that normally available from the front panel. These 
functions are intended for the user who knows the instru- 
ment and the service technician who needs arbitrary con- 
trol of the instrument functions, Automatic tuning and 
ranging, over voltage protection, and error messages protect 
the user from invalid measurements during norma] opera- 
tion. When special functions are used, some of these 
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safeguards are removed, depending on ihe special function 
selected, and thus there is a degree of risk that the mea- 
surement may be invalid. However, there is no risk of dam- 
age to the instrument. 

To enter a special function* the user enters the special 
function code (usually a prefix, decimal, &nd suffix) then 
presses the SPCL key. The special function code appears on 
the display as it is being entered. If a mistake is made during 
entry of the special function code, the user can press the 
CLEAR key and start over. When a special function is en- 
tered, the light in the SPCL key goes on if it is not already on. 
The readout on the display depends on the special function 
entered. It may be a measured quantity, an instrument set- 
ting, or a special code. In some cases the display is unal- 
tered Special hi net ions can be entered from the HP-IH by 
issuing the special function code followed by the code SP, 

Floating Input and Output 

To eliminate troublesome ground loops, both the source 
output and the measurement input of the 8903 A are \l -i- 
ing, This is helpful in tow- distortion or low-level ac mea- 
surements when it is necessary to reject potential differ- 
ences between the chassis of the B9Q3A and the device 
under test. The 8903 A also has EMI [ electromagnetic, inter- 
ference] protection built into the source output and mea- 
sure men I input lines so that it can work in the presence of 
high EML All of the analog circuitry is shielded by an 
internal LAll-tight box. The output and input lines extend- 
ing from this box to the front panel are shielded and termi- 
nated in BNC connectors. For user convenience, the BNC 
connectors are spaced so that BNC-tn-banana adapters can 
be attached. Thus a banana or twisted- wire connection can be 
made to the instrument when EMI shielding is nol critical, 

Analyzer Architecture 

The H 903 A Audio Analyzer combines three instruments 
into one: a low-distortion audio source, a general-purpose 
voltmeter with a tunable notch filter at the input, and a 
frequency counter. Measurements are managed by Ihe 
microprocessor-based controller. This combination can 
make most common measurements on audio circuits au- 
tomatically. To add to its versatility, the analyzer also has 
selectable input filters, logarithmic frequency sweep, X and 
Y outputs for plotting measurement results versus fre- 
quency, and HP-IB programmability. Fig. 2 is a simplified 
block diagram. 

The amplitude measurement path flows from the INPt'T 
jacks [HIGH and LOW) to the MONITOR output on thi 
panel, and includes the input and output tms detectors, the 
dc voltmeter (the voltage-to-time converter and counter), 
and the SIN AD meter circuitry. Measurements are made on 
the difference between the signals at the hk;h jaok and the 
LOW Jack. Differential levels can be as high as 30CJV. Signals 
that are common to both the HIGH and LOW jjick.s are bal- 
anced out. Signals applied to the LOW jack must not e\- 
i eed 4V 

The input signal isac coupled for all measurement modes 
except dc level The signal is scaled by the input attenuator 
to a level of 3V or less. In pruteei the ,u live circuits, the 
overvoltage protection circuit quickly disconnects the 
Input amplifier if its input exceeds 15V 



The differential signal is converted to a single-ended 
signal (referenced to ground) and amplified. The sig- 
fiuiher amplified by a programmable gain amplifier, which 
is ac coupled. The gain of this amplifier and the 
differentia I -to-single-ended amplifier are set to keep the 
signal level at the input rms detector between 1*7 and 3V 
rms to optimize its effectiveness and accuracy 

The output from the first programmable gain amplifier is 
converted to dc by the input rms detector and measured by 
the dc voltmeter. The output of the detector is used to set the 
gain of the input circuits and becomes the numerator of the 
SIN AD measurement and the denominator of the distortion 
measurement. The gain of the input path is determined by 
measuring ihe dc level. The input rms detector also 
monitors the ac component (if there is one] and lowers the 
gain of the input path if the ac signal will overload the input 
amplifier. At this point either the 400-Hz high-pass filter or 
the psopho metric filter can be inserted into the signal path. 
The 400-Hz high-pass filter is often used to suppress line 
hum or the low-frequency squelch tone used in some 
mobile receivers. The psophometrie filter has a bandpass 
frequency response that simulates the "average" response of 
human hearing, It is often used to condition a receiver audio 
output when determining the receiver's Input sensitivity. 

During SIN AH. distortion, or distortion level measure- 
ments, the fundamental of the signal is removed by the 
notch filter The output from the filter is the distortion and 
noise of the signal. In the ac level and signal-to-noise 
modes the notch filter is bypassed. After amplifying and 
low-pass filtering, the output from the notch filter is con- 
verted to dc by the output rms detector and measured by 
the dc voltmeter. 

During distortion or distortion level measurements, the 
notch filter is tuned to the frequency counted at its input. 
Coarse tuning is done by the controller, and internal analog 
circuitry fine tunes and balances the notch filter. During 
SINAI) measurements, the controller coarse tunes the notch 
tn the source frequency, Thus a SINAI) measurement is 
normally made with the internal source as I he si i mill us; this 
permits measurements in the presence of large amounts of 
noise (where the < < mil oiler would be unable to determine 
the input frequency). If an external source is used in the 
SINAI) measurement mode, the source frequency must 
be within 5% of the frequency of the internal soup e 

The two programmable gall) amplifiers following the 
notch filter amplify the low-level noise and distortion sig- 
nals from the notch filter. The overall gain of the two 
am p\ i Tiers is normally set to maintain a signal level of 0.25 
to 3V at the output detector and monitor. The ;4()-kHz ami 
80-kJfz low-pass filters are selected from the keyboard. 
With nn low-pass filtering, the bandwidth of Hie measure- 
merit system is 7 50 kHz, The filters are most often used tn 
remove the high- frequency noise components in low- 
frequency distortion and signal -tn-m use measurements. 
The output from the second programmable gain amplifier 
drives the rear-panel MONITOR on I put jack. Taking advan- 
tage of the increased amplification available at this point, 
ounter monitors this output in ac level and signal-ln- 
00136 modi 

The output rms detector is read by the dc voltmeter in the 
ac level. SINAI) (the denominator), distortion [the 
numerator), distortion level, and signal-tCHnoise measure- 
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Audio Analyzer Applications 



The 8903A Audro Analyzer's measurement capabilities reach far 
beyond conventional distortion analyzers. Much of this performance 
results from microprocessor control and HP- IB programmabiJity 
Numerous hardware features such as a fast counter, both analog 
meter and digital display, and switchable detector filtering allow the 
user to make unusual or special measurements with convenience 
and little auxiliary apparatus 
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Consider the 8903A used at a test and calibration station in the 
manufacture of audio power amplifiers A typical sequence of events 

mighi include an output offset null, frequency response check, distor- 
tion test, and noise measurement , The 8903A can perform all These 
measurements qupckty with a single test setup If an X~Y plotter is 
connected to the rear-panel outputs, the results of swept frequency 
measurements can be recorded on standard log/lag or log/lin 
graph paper. 
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Fig. 1. Stvepf distortion and fre- 
quency response of a two-pole ac- 
tive fitter, measured by the 8903 A 
Audio Analyzer. 



ment modes. It is also used to set the gain of the two pro- 
grammable gain amplifiers- The (Intent or can be configured 
internally to respond to the average absolute value of the 
signal instead of the true rms value. This option is provided 
because some measurement specifications for detection of 
distortion and noise specify the use of an average respond- 
ing detector, Average responding detectors do not give an 
accurate indication ut signal power unless llu- si^n.il wave- 
form is known, I If the waveform is Gaussian noise the read- 
ing will be approximately 1 tlB less titan the true rms value.) 

In the SINAD mode the outputs from the input and output 
rms detectors are converted to logarithms, subtra< ted, and 
converted to a current by the SINAD meter amplifier to drive 
die SINAD panel meter. Since SINAI) measurements are 
often madr under very imisy conditions, the panel meter 
makes it easier to average the reading and to discern trends, 

The voltage-to-time converter converts the dc inputs into a 
time interval, which is measured by die counter. 

The 89-03 A uses a reciprocal counter. To measure fre- 
quency, it counts the period of one or more cycles of the 
signal at its input, Then the controller divides the number of 
cycles by the accumulated count. The reference tor the 
counter is the 2-Vihbc time base, which also is the clock for I he 



controller. The counter has four inputs and three modes of 
operation: 

1. Voltage measurement. The time interval from the 
voltage-to-time converter is counted. The accumulated 
count is proportional to the dc voltage, For direct mea- 
surements (ac level and distortion level], the count is 
processed directly by the controller and the result is dis- 
played. For ratio measurements [SINAD, distortion, and 
signal-to- noise), the counts of two successive measure- 
ments are processed and displayed. For SINAD and dis- 
tortion, the controller computes the ratio of the outputs of 
the input and output rms detectors. For signal-to-noise 
measurements, the output of the output rms detector is 
measured with the oscillator on and off and the ratio of the 
two measurements is computed. 

2. Input frequency measurement, The signal from the last 
programmable gain amplifier or the high- pass bandpass 
tillers is conditioned by the counter input Schmitt trig- 
ger to make it compatible with the counter's input. The 
period of the signal is then counted, the count is pro- 

i d by the controller, and the frequency is displayed, 

3. Source frequency measurement. The counter measures 
the frequency of the oscillator during tuning and when 
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Fig. 2. est etup for screening operational amplifiers 

The swept frequency measurement with plot capability finds many 
applications in the laboratory Fig i shows the swept distortion and 
frequency response of a two-pole active filter The upper curve shows 
the magnitude response, while the tower shows distortion. Notice that 
the analyzer input magnitude covers a 30-dB dynamic range. During 
the sweep, the 8903A is automatically setting input gain before per- 
forming the distortion measurement. 

If an H P- IB controller is available, mere are many more applications 
for the 8903 A. Fig. 2 shows a simple test setup for screening opera- 
tional amplifiers With no other instruments in the system, a 
computer- controlled 8903 A can rapidly and accurately measure 
input offset voliage. input noise voltage, and distortion It can also be 
used to measure the gam- bandwidth product of the op-amp, pro- 
vided it is nol greater than 30 MHz. The controller test program can be 
written to provide either a go/no-go output or a listing of measurement 
res u Its - 

Many of the special measurement modes are available through the 
use of special functions. For example, the 8903A can be used as a 
test amplifier with gam settable from -24 dB to +94 dB. Signal 
filtering can be added by selecting the appropriate front-panei con- 
trols, or the special functions can be used to put the instrument into a 
notch or bandpass filter mode. Of course, signal frequency and 
amplitude will be measured and displayed if the operating loads are 
chosen correctly, Th*s mode of operation makes it possible to count 
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Fig, 3. Transmitter test setup ustng the 8901 A Modulation 
Analyzer and the 8903 A Audio Analyzer 

the frequency of a signal whose amplitude is in the low tens of 

microvolts. 

Transceiver Testing 

Much effort has gone into the design of the 8903A to facility I 
audio measurements required for automatic, programmed. 
transceiver testing. For example, the worst-case source frequency 
error of 0.3% allows the 8903A source to replace a synthesizer for 
squefch-tone generation In addition, a binary programming m 
available th ro ug h the HP- 1 B that causes the 8903 A to generate a tone 
burst sequence that can be used to unlock a coded receiver A 
related function allows the 8903A to measure burst tones generated 
by an external source, such as a transmitter under test. The packed, 
four-byte output allows the 8903A to output Frequency measurements 
as often as every eight milliseconds. 

When the 8903 A is used in conjunction with the 8901 A Modulation 
Analyzer, almost all transmitter tests can be automated. Fig. 3 shows 
the block diagram, With the source turned off and the transmitter 
keyed, the squelch-tone frequency can be counted Then the 
400-Hz high-pass filter can be switched in to eliminate the squelcn 
tone from the remaining measurements. With the source output turned 
on. the 8903A can easily be programmed to make the necessary 
measurements to determine distortion and microphone sensiti. 



verifying that the oscillator frequency is within toler- 
ance. This frequency is normally not displayed 
The source covers the frequency range of 20 H/ In 100 
kHz, It is tuned by the controller to the frequency Bill 
from the keyboard, using a tune-and-couM rout i he:. (Mole 
that the frequency is not obtained by frequency synthesis.) 
The switch following the oscillator is closed except in the 
signal-to-noise measurement mode or when an amplitude 
of 0V is entered from the keyboard. 

The source output amplifier and output attenuators pro- 
vide 77.5 dB of attenuation in 2,5-dB steps, and 2.5 dB of 
attenuation in 256 steps. This gives an open-circuit output 
from 0-0 niV to hV\ The filiating output amplifier converts 
the ground-referenced input to a floating output. Either 
output, HIGH or LOW, can be floated up to ten volts peak. 
The entire operation of the instrument is under control of 
the microprocessor-based I cinl roller, which sets up the in- 
strument at turn-on, interprets keyboard entries, executes 
changes in the mode of Operation, continual I v monitors 
instrument operation, sends measurement results and er- 
rors to the front-panel displays, and interfaces with the 
f IP-IB, Its computing capability is also used to simplify 
circuit operation. For example, it forms the last stage of the 



counter, converts measurement results into ratios [in % nr 
dH|, and soon. It also executes rmi lines for servicing the rest 
nf the instrument as well as its* -II 

Input Circuits 

Numerous design constraints were fmposed on the input 
attenuator/protectionamplilier circuitry. An input imped- 
ance of 1 00 kO is necessary to prevent the input circuitry 
from unnecessarily loading the device under test. On the 
other hand, maintaining a good signal-to-noise and distor- 
■i on ratio, good frequency response to ion kHz, automatic 
Operation and reliable performance with input signals as 
large as 300V is very challenging. Consider the 300V con- 
st r,nnt. It is possible for a 300V signal to appear suddenly at 
the input while the instrument is measuring a 50- mV input 
Level. Mot only must no damage occur* but also the overload 
recovery must be quick H and the input protection circuit 
must not be all n wed to degrade I he input noise floor when 
measuring 50 mV. The 300 V signal may also have spikes or 
transients that lar exceed 300V, or the user may inadver- 
tently apply a larger signal- in all those cases the circuit 
must recover without causing a safety hazard to the user, 
destroying internal components, or even blowing a fuse. 
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Making the Most of a Microprocessor-Based 
Instrument Controller 

by Corydon J. Boyan 



In the 8903A Audio Analyser, most of the tasks that lead to the 

display of measurements are coordinated by a microprocessor The 
processor {an F8 with 18K bytes of ROM and 192 bytes of RAM) 
counts and tunes the interna! audio source, sets the input amplifier's 
gam. tunes the notch filter, sets the output amplifier's gain, and 
measures the voltages that will be used to generate each reading 
This means that the B903A can, among other things, automatically 
take distortion readings without calling upon the user to turn several 
knobs when seeking a null Th<s ability alone is ample justification for 
basing the instrument controller on a microprocessor, but the 8903A 
goes tar beyond this in applying the power of its controller 

"Guaranteed' 1 Accurate Measurements 

The 8903A is HP- IB programmable, and this brought some impor- 
tant factors into consideration during the design. The performance of 
the internal source when settling from one frequency or level to 
another, for example, is a key factor in assuring the validity of the first 
measurement taken after changing the frequency or level. Similarly, 
the settling performance of the input and output amplifiers and the 
notch filter under various adverse signal conditions becomes very 
important when the instrument is trying to deliver an accurate first 
measurement after a change in operating parameters (e.g., after 
tuning to a new frequency). Add to this the desire to make measure- 
ments as rapidly as possible and you have a very interesting problem 
for a microprocessor to soive, 

For example, every lime the internal source is tuned, the processor 
spends about 1 70 mHltseconds tuning to and verifying the frequency 
This translates to over 66,000 operations. Of equal complexity is the 
pb of setting the correct gam for the output amplifier and allowing the 
circuits to settle before the output amplifier voltage is read, with the 
object of never giving the user an invalid reading while at the same 
time delivering the reading in as short a time as possible. To ac- 
complish this, the routine that controls the output amplifier makes use 
of such techniques as slope checking and frequency dependent 
delays to ensure rapid, valid readings. These techniques make pos- 
sible accurate readings in half the time it might otherwise take. 
Referring to the flow chart, Fig 1, note that the key to this routine is the 
technique of measuring the rate of change of the signal on the output 
amplifier (the slope) before attempting to check the signal level to 
determine if the output amplifier gain is properly set This is done 
because it is quite common, after the signal level has suddenly 
changed (for example, when the notch filter is suddenly switched in 
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Fig. 1 , 8903 A leveling algorithm. The microprocessor checks 
the rate of change of the input signal as weft as its tevet to make 
certain that the signal has settled before the output amplifier 
gam is set. 

when going from ac level to distortion mode), to have a rapidly falling 
output amplifier voltage pass through the acceptable region, and 
thus have the gain appear to be properly set. when in fact the signal is 
on its way to a level that wilt require more gam for the proper reading 
To keep the controller from being fooled by this phenomenon, we take 
two voltage readings in rapid succession, and from them calculate 
the rate of change of the signal. The time over which we measure this 
slope is longer at lower frequencies, and is in fact looked up in a table 
in ROM, based on ihe frequency at which we're operating II the rate 
of change is too fast during this period, we delay before checking the 
level (see Fig 2). This delay is also frequency dependent and table- 
driven. Note that if we simpty delayed each cycle to keep Irom risking 



Dc inputs and low-level differential amplification for 
comnion-mode rejection are necessary features for an in- 
strument likn the 8903 A. One consequence of automatic 
ranging is that low-capacitance mechanical switching 
techniques cannot he used effectively. Needed are high- 
voltage reed relays, which affect high- frequency perfor- 
mance and require the use of compensation capacitors, 

For dc operation the first part of the input signal path isdc 
coupled with the input blocking capacitor bypassed. The 
output of the different ial-ta-single-endfid amplifier can 
then be monitored to obtain an accurately scaled represen- 
tation of the input dc level 



input voltages larger than 3V are attenuated by the input 
attenuator, a network of resistors that divide down the input 
signal. The appropriate tap point ts selected by a reed relay 
network, If an overload occurs, the maximum attenuation 
setting is enabled. 

To protect the sensitive input amplifier following the 
attenuator from short-term transients, an overvoltage pro- 
tection network is used. For low-level signals the transfer 
impedance is low and signals applied to the input are 
coupled to the differential-input amplifier. However, for 
input signals large enough to damage the amplifier the 
output of the protection network is limited to a safe value. 
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Signal while settling 



■ Level looks good here 



k 
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v 
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but slope measured in this interval 
shows need lo delay to let signal settle 



Fig. 2. An example of slope checking by the 89Q3A micro- 

processor 

displaying invalid readings, every measurement would take longer, 
perhaps by several hundred milliseconds. The slope checking gives 
us a rapid check on signal quality which eliminates the need to delay 
every June 

After the slope check (and possible delay) we assume the output 
amplifier is well settled and ready for leveling. We then measure the 
voltage at its output, and if too high or low. adjust the gain accord- 
ingly. After adjusting the gain, we make another pass through the 
leveling algorithm (unless we are now \n the highest-gam range). To 
keep from getting caught in infinite loops, we require that after the first 
pass the gain never be reduced, and assume that if it is, we have an 
unstable signal, which causes us to start the entire measurement 
cycle (tune source, level input, tune notch, level output) over again 
The final protection from infinite loops comes from counting the 
number ot times we restart without displaying (each time we put a 
"— -" pattern on the dispfay). ana after 128 times we display Error 31 , 
which also goes out to the HP-IB 
Sweep 

Because the 8903A has the ability to generate accurate and com- 
paratively rapid measurements automatically, with no need for the 
user to insert delays m HP-IB routines or wail for the display to settle 
before taking a reading, it is possible to have the instrument sweep 
itselt over a range of frequencies without user or computer interven- 
tion. Here the calculating power of the microprocessor is brought to 
bear on the problem of determining the frequency increment for each 
new point in the logarithms sweep, based on the sweep range and 
the number of points per sweep (which can be set by the user). All the 
user need do to get a series of measurements spanning a frequency 
range is set the start and stop sweep frequencies and the number of 



points per sweep, and press the nakes it easy to 

jre the frequency response of an amplifier 

X-Y Recorder Output 

The calculating power of the microprocessor-based controller s 
also apparent in the operation of the 8903A s X-Y recorder outputs 
These e driven by dsgrtaMchanalog converters controlled 

bytr>emicroprocessor.ratherthand rs As 

a result, the user does not have to worry about the recorder output 
voltage abruptly changing when the analyzer autorange 
era processor scales the recorder output according to the displayed 
reading and the plot limits entered via the keyboard. The recorder 
outputs are always between zero and ten volts so the recorder s zero 
and vernier controls need be adjusted only once Thus a properly 
scaled plot is easjly generated by using the sweep and the X-Y 
outputs, without any need for an external controller 

Special Functions 

Hidden behind the basic measurements are nearly forty special 
functions, whrch provide extended measurement capability and 
many service aids. For example, the analyzer can be given a load 
resistance in ohms and commanded to display ac level m watts 
Another special function changes the number of points per decade in 
a sweep, and several special functions modify display operation. 
Normally the left and right displays indicate the frequency and level 
(or distortion, etc ) of the srgnal applied to the analyzer input. Some- 
times, when using the 8903A just as a source, the user may want the 
analyzer to display the frequency and level of the source Special 
function to provides this display 

Service aids provide front-panel display of many internal voltages 
and settings. Without microprocessor control each special function 
would require one or more switches on the front panel instead of one 
special key, and would therefore probably not be included Thus, the 
processor allows implementation ot useful features the user would 
not otherwise get. 
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The network consists of two back -to- back diodes which 
open up under large-signal conditions. 

The differential amplifier consists of three high- 
performs n< t 1 amplifiers These have the necessary noise* dr. 
offset and frequency performance so that they do not de- 
grade the signal quality. The total input amplifier chain acts 
as a 4-dB/step amplifier with leveling hysteresis, which 
maintains the post-amplifier signal level within 6 dti (3 to 
1.5V rms). Should the output level change with time and 
deviate from this range, the gains of the attenuator, the 
diffcrirntial-tn- single-ended amplifier, and the program- 
mable gain amplifier are adjusted to compensate. 



To summarize, the gain of the input amplifiers is mod- 
ified in three ways, First, the microprocessor monitors the 
input detector. If the detector voltage is too high or too low. 
the microprocessor varies the gain of the attenuator 1 ' 
amplifier chain to bring the level vvilhin bounds. Second, 
the overvoltage protection network limits if the input signal 
exceeds ±15V, Third, any sustained input overload trips 
the input overload detector- This detector monitors the 
input rms detector and the differential amplifier, and if 
either exceeds a certain absolute voltage limit, the overload 
detector trips, resetting the gain of Hn> > niin j amplifier 
chain to its minimum value (maximum attenuation). This 
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Design for a Low-Distortion, Fast-Settling Source 

by George D. Pontis 



To fulfill, the requirements of the 8903 A Audio Analyzer ihe built-in 

source must have good performance in certarn key areas. First of all. 

for swept measurements, it must be readily programmable and fast 

settling. It must ateo have Sow distortion and noise, and it must have 

jbod arnpJrtude accuracy over ihe entire frequency range of 20 

100 kHz. 

The combination ol these conflicting requirements suggested the 
use of a high-performance RC oscillator instead of a synthesizer The 
synthesizer is easy to program and settles quickly, but it >s difficult to 
bund a synthesizer with noise and distortion more than SO dB oelow 
the fundamental Synthesized designs also do no\ nave sufficient 
absolute level accuracy or flatness without leveling, and a leveling 
loop that does not unduly degrade the distortion and settling time 
would be very difficult to desjgn. 

Unfortunately, none of the common RC audio oscillator cfe&grfis 
looked suitable either Usually tne amount of distortion is inv- 
proportional to the settling time. Also, the tuning elements usually 
float making tne circuit dillicultto interface with programming lines. 

For these reasons a state-variable oscillator was cnosen. Sfrnilar to 
that proposed by Smith and Vannerson in 1975 1 Since this oscillator 
is built around a state-variable filter structure, inexpensive JFET 
switches can be used easily to switch the tuning elements More 
important is that the ALC design permits very rapid settling without 
tracing oft good distortion performance. Although the programmed 
frequency does not have the accuracy of a synthesizer sufficient 
resolution is available to permit firmware tuning to within ±0.3% of 
the programmed value 

Tnis oscillator design can be described as a state-variable filter in 



h the resistor that determines the Q is replaced by an analog 
multiplier Tne control signal for the multiplier is provided by an 
automatic level control {ALC) circuit. The ALC circuit compares the 
oscillator amplitude with a stable dc reference obtained from a 
temperature-compensated reference diode The resulting error sig- 
nal is processed tnrougn two paths One path carries tne cycie-oy 
cycle proportional error to the controlling multiplier This provides 
very last setting The other path includes an integrator in the loop to 
eliminate nearly all of the steady-state error Tnis design is theoreti- 
cally capable of settling the output amplitude within two cycles after 
small-signal amplitude disturbances 

There are two important refinements in the 8903A oscillator The 
first is the use of a special two-stage peak detector This consists of 
track-hold and sample/hold amplifiers to eliminate any distortion- 
causing ripple on the detected peak output The second refinement is 
the addition of an ALC loop gain control to compensate the leveling 
loop, cycle by cycle, for changes in oscillator amplitude This greatly 
decreases the large-signal settling time 6f the oscillator, which s 
important when switching from one range to another 

Fig 1 shows the oscillator integrators. The gain constant of these 
integrators is changed in three-octave steps by selecting the feed- 
back capacitor. This gives us range switching. Coarse tuning within 
each range is done by switching the input resistors. The eight 
binary-weighted resistors provide 255 usable steps Psacing the 
switching devices at the virtual gound point permits the use of JFETs 
with 3ow dram -to- source on resistance (R D I An individual transis- 
tor switch conducts when its gate is allowed to rise to ground poten- 
tial, and turns oft when its gate is pulled to the negative supply, -15V 
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Fig. 1, The gain-switched inte- 
grator used in the 8903 A' s inter- 
nal oscillator Capacitors are 
switched to change ranges Resis- 
tor switching provides coarse tun- 
ing within each range 



n^ot occurs within ten milliseconds ii ihe overload i 

This protects the input from burning out or blowing a 
ttid allows for rapid OV& Icrad 



Notch Fitter 

The notch filter design challenge was twofold, First, it 
was necessary to design a lmv-distortion< low-noise filter 

that was also programmable. Second, to minimize overall 
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The switching t . 

-3es and quae "- .- to TTL 

levels. 
Fig 2 is a b'C itor Integrators Ui and U2 

and jnven- 

done by (J -- ' hed resistor nc 

; 
switched tne gam of LM changes, effectively altering the am:. 
signal transmitted frc The gam is proportu 

\ --BH, where R ' e selected input 

resistors and A and B are constants fine tuning 

range. This gives the oscillator enough resolution to tune 
=t 2% of any frequency within the range of The instrument. 

11 t s generally true for sinusoidal oscillators that purity and settling 
time are chiefly functions ot the ALC circuits or mechanisms used 
The two-stage peak -detector circuits are the key to the performance 
of the 8903A oscillator Oscillator amplitude data is obtained by the 
f • - hold and sample hold amplifiers in the following manner Switch 
S1 closes during the time the output is at its negative peak. Capacitor 
C1 rapidly charges following the sine-wave amplitude up to its posi- 
tive peak At thus time S1 opens, holding the peak voltage on CI 
Switch S2 then quickfy closes and opens again, thus updating the 
sampled peak level held by capacitor C2 

The two-stage scheme has several advantages For one, the first 
stage may be optimized for fast da! a acquisition, while the second 
state is optimized for long hold lime. Fast acqi itfal for 



gooc : - _ '. : 

lam r ^wd«sto r " so thissel 

nas no stea 

■ appear output a r mam 

■ 
■ 

functions are accom: 
gam buffer 

Fo f leral. settling t»me is stn 

dependent upon loop gair . 

r proportional to f re : : oy switch- 

er integrator resistor on once per cycle for a duration of 35 ^s 
'eases the integral error signal with frequency The 
duty cycle increases with frequency until the integrator switch S3 is 
left closed continuously for frequencies greater than 25 kHz Above 
25 kHz, the oscilfator easily settles in less than one millisecond 

3903A osciflator performance is largely limited by the quality of 
commercially available, reasonably priced analog multipliers The 
multiplier used is decoupled slightly (from optimum) to reduce its 
contribution to THD and noise This extends the oscillator settling time 
to a period of four to five cycles 

Reference 

■ ■ 

■ 




4,24V 
Oscillator q V 

Output 

-424V 



Fig. 2. Block din , ' -e 8QQ3A\ mt&rfial state-variable oscillator 



measurement time, it was necessary to develop an i 
rapid tint;- tune mechanism with quick recovery from over- 
loads and mistiming. To achieve good Q over the frequeni $ 
range, an active RC filter ts necessary. Ta tune this device, a 



variable resistance ar conductance device ts needed, since 
achieving the tuning range with variable capacitors ur in- 
ductors is tapractii al. Many resistiveiy tunable a< live < on- 

figurations are feasible- However, determining the up- 
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Floating a Source Output 

by George D. Pontis 



To provide the greatest versatility m both benchtop and systems 
applications, the 89D3A Audio Analyzers built-in source *s floating 
This lets the user eliminate ground-loop errors, sum signals, and add 
dc offsets to the source output 

Previous designs used a separate, isolated power supply for the 
source circuits. This method is straightforward and offers very good 
Low-frequency common-mode rejection, However there are several 
reasons why this arrangement is not used in the 89Q3A 

The biggest problem with the floating power supply approach is 
interfacing with the digital programming lines Since only three of the 
thirteen attenuator lines have relay isolation and none of the nineteen 
oscillator lines can be floated, over thirty lines must be coupled in 
some manner. One solution is to float on«y the final output stage. This 
eliminates the need tor couplers, but requires a high-performance 
differential-input amplifier to reject the common-mode signal that 
appears at the inout of the floating stage. Since a floating power 
supply is still required, the cost of this approach is relatively high 

The 8903A solves this problem with a single-ended-to-differential 
output converter This circuit shown in Fig. 1 , operates on the instru- 
ment s ground-referenced ±15V supplies and requires only two op- 
erational amplifiers. A precise combination of negaiive feedback, 
positive feedback, and cross-coupling yiesds a symmetrical differen- 
tial output vvith infinite common- mode rejection and a well-defined 
output impedance 

An analysis of this circuit is generally a tedious procedure because 
of the number of components involved. However , the high degree of 
symmetry m the circuit can be exploited to great advantage by using 
the relations R2/R1 = R12/R11, R3 = R7 R6 = RIO, and R4 - R5 = 
R8 = R9. From these relationships, one can derive the expression 
R2/R1 = (2R6 + R4)/2R3, which is a necessary condition for achiev- 
ing an infinite common-mode output impedance. Then it is easy to 
calculate the differential output impedance and the open-circuit volt- 
age ga«n. The resulting equations can be manipulated to find suitable 
values. For the resistor values used in the 8903 A, the associated 
gam is 1 125 and the output impedance is 480 ohms. The output 
is further padded with a 120-ohm resistor to yield the desired 600- ohm 
output impedance. 

It would have been possible to use resistors that gave an output 
impedance of exactly 600 ohms instead of 480 ohms, but this would 
have required setting up and stocking a supply of several extra odd 
resistor values. As it is, the circuit is realized using Q.1%. 25-ppm 
resistors that are also used elsewhere in the instrument 



tirnum match between the filter configuration and the vari- 
able resistive element is not straightforward. 

Let's go through the alternatives and the tradeoffs, 
Switchable resistor networks have good distortion and 
noise characteristics but do not provide continuous tuning 
coverage and require extensive switching circuitry Photo- 
resistors can he driven over a large resistance range and 
provide continous tuning However, the noise and distor- 
tion they add to a signal are greater than the required level 
of 90 dB below the signal level. They can be used as fine- 
tuning elements if coupled only partially into the circuit, 
These devices can also be slow and are awkward to control 
mp idly* reducing the tuning speed. Finally, they tend to 
vary significantly from device to device and with time and 
temperature, making compensation difficult. 
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Fig. 1 , Singie-ended-to-differentia! output converter provides 
a floating output for the 8903 As internal Qsdliatci 

The easiest way to see how this circuit works is to eliminate either 
the inverting (low) or noninverting (high) half of the circuit by shorting 
the respective output to ground. Fig 2 shows the reduced circuit 
when rhe low half is grounded, If R4 is disconnected, the circuit will 
have a forward gain of about one. and an output impedance oi 276 
ohms. R4 works in conjunction with R6 to provide voltage and current 
feedback that causes the gain and the output impedance to rise. 

To demonstrate that the output js truly floating, we ground the input 
and apply a test source to both outputs. Ideally, the current flow from 
the test source should be zero. Fjg. 3 shows a block dsagram and the 
reduced circuit for this test Here it can be quickly calculated that tne 
output of U1 will rise just enough over that of the test source to make 
the current through R6 cancel the current through R4 and R8. Note 
that the current flowing through sources V1 and V2 is supplied by the 
other halt of the circuit, which is not shown. 



Four-quadrant analog multipliers also do not have the 
90-dB performance necessary, but they too can be lightly 
coupled into the circuit for fine tuning. These devices are 
fast, inexpensive, and easy to drive, There are many varia- 
tions on this type of circuit, some of which can he obtained 
in. integrated form, Those most suitable use a differential 
pair of bipolar transistors as a variable gain element by 
varying the common -mode current, 

Light bulbs as variable resistive devices are relatively 
linear but are slow and have a limited, dynamic range, 
Thermistors, diodes ami other nonlinear devices would all 
be useful only for fine-tune applications. The drive and 
compensation circuitry for all of these alternatives would 
be complex and the overall performance marginal. 

The tuning elements selected were swi tenable resistor 
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Fig. 2. Eliminating the inverting side of (he circuit of Fig 1 Dy 
shorting the low output to ground results m this reduced cir- 
cuit. 

In practice it was found that parasitic effects and the 
electromagnetic-interference (EMI) filters degrade the circuit bal- 
ance when the frequency approaches 100 kHz, However, each 
board is tested for a minimum of 50 dB common-mode rejection at 1 
kHz A typical unit has greater than 40 dB refection at 100 kHz. Ai so, a 
1 0-fcO resistor internally i les t he tow output to the chassis g round . This 
provides a reference for the output when no external load is con- 
nected. Without Ehis resistor, the common-mode output voltage is 
indeterminate. 

One initial concern about the circuit was difficulty of troubleshoot- 
ing problems, such as one of the resistors drifting out of tolerance, 
causing poor common-mode rejection. This problem was solved by 
implementing the following test procedure. First, the input to the 
Moating ampli-ier is set to exactly 1 00V mns using the special func- 
tions built into the 8903A. Then, the technician shorts the low side to 
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Fig. 3, fa j To demonstrate that the output of the circuit is 
floating, the input is grounded and a voltage source is applied 
to both outputs . (b) Reduced circuit for this test (only half of the 
arcutt is shown). Current through R6 cancels the current 
through R4 and RB Current through V1 andV2 is supplied by 
the Other half of the circuit (not shown) Thus no current is 
drawn from the test source 

ground and measures the potentials al several circuit nodes 1 1 e 
measured values can be compared to the calculated values pub- 
lished m a table in the S903A manual A second set of measurements 

can be made, if necessary, with the high side grounded, By observ- 
ing the deviations between the measured and calculated values, it s 
easy to locate the faulty component. 



networks for low-distortion coarse tuning, and a four- 
quadrant multiplier for fine tuning between the discrete 
steps of the swif enable resistors. Since the four-quadrant 
multiplier is coupled into the circuit only enough fur ±7% 
tuning, it does not contribute significantly to the overall 
noise and distortion. A resistor switching network is best 
implemented if one end of the network is dynamically kept 
at ground potential; this relieves many constraints on the 
switching network. To this end, a state- variable* notch 
configuration is used (see page 14). With this design, low- 
distortion tuning over a 10:1 frequency range is achieved 

tf&to'Variablti rUftf IS a consequence ol the lacl that the equations thai de 
se n be ih-e system cari be ,v ft! ,!uj ... upm:' . 

aJnearayetern, Lis * Ax Bu - « liable) oMhs system 

lo an input u 



Capacitors are switched into the network to change fre- 
quency in three- octave bands and provide complete cover- 
age ol the frequency range of the analyzer. 

To complete the fine-tuning path , a synchronous detector 
mixes the filter output waveform with the fundamental 
waveform. If any fundamental exists in the notch output, a 
dc: current is generated to fine tune the notch. The critical 
parameters here are 100-dB dynamic range and rapid opera- 
tion. A PET double-balanced mixer was selected. The input 
mixing signal that drives the FET is a square wave, rich in 
odd harmonics, so the circuit responds to odd-order har- 
monics as well as to the fundamental, This is the classical 
solution. A complete null may not be achieved if third, fifth, 
or higher-order odd harmonics wte present. Total error can 
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A Digitally Tuned Notch Filter 

by Chung Y, Lau 



The notch filter in the 8903 A Audio Anaryzer rejects the fundamen- 
tal frequency component of the incomrng signal. This fitter consists of 
a state-variable active filter and fine-tune and fine-balance control 
circuits. A simplified schematic of the filter minus the control circuitry 
is shown in Fig. 1. 

Conventional notch filters (bndged-T, etc.} used in distortion 
analyzers are not well suited for digital control because they require 
expensive relays or analog switches. The state-variable filter ap- 
proach has the following advantages 

The Q of the filter is hxed and independent of frequency. 
Tuning ts accomplished by switching resistors and capacitors in 
the two integrators. Inexpensive JFETs are used as switching 
elements and the switch drivers are simple because both the gate 
and the channel of the JFET are virtually at ground potential when 
the switch is on. 

Three filter outputs (low-pass, bandpass, and high-pass) are 
available simultaneously. The low-pass and high-pass outputs are 
phase-shifted from the input by 90" when the filter is tuned to the 
input frequency. These two outputs are used in the control circuits, 
so extra phase-shifting circuits are not needed 
Distortion generated in the fine-tune circuitry is filtered by two 
integrators before appearing at the notch output 
In distortion mode, the filter is tuned in the following manner The 
microprocessor counts the input frequency and tunes the notch filter 
to the same frequency by switching In the proper capacitors and 



resistors When the filter is tuned to the fundamental frequency of the 
input signal the bandpass output {\fa ? in Fig, T) inverts the funda- 
mental component of the input signal U5 sums V, N and V B: ; to cancel 
the fundamental component of V irst exactly However, because of the 
phase and amplitude characteristics of V BP , the harmonics present in 
V r , are relatively unattenuated. 
Analytically, we have 



— (s) 



-s*jjQ 



S L + SmJO 4 tiif, 



and 



V np 4V lh 



'(*) = 



S 3 4 SWw/Q + (ill 



where a» Q is the center frequency in rad.'s, Therefore, there is zero 
transmission at s = ±j^,, 

The state- variable filter alone can provide only about 15 dB of 
fundamental rejection because of the coarseness of the frequency 
tuning (use of discrete resistor values, resistor mismatches, etc.). 
Fine-tune and fine-balance control circuits are necessary to achieve 
a notch depth greater than 90 dB. The fine-tune circuit insures that the 
filter is exactly tuned to the input frequency (no phase error) and the 



Correction Signals 
from Control Circuits 



Input 




Fig. 1. 8QQ3A notch filter. Correc- 
tion signals from the control cir- 
cuits of Fig. 2 provide fine tuning 
and fine balance. 



be as much as G\46 dB. but this is considered reasonable. 
Rapid response is achieved in two ways. First, when an 
overload occurs, the circuit gain is quickly reduced by the 
output amplifier circuit, minimizing the transient and thus 
minimizing the recovery time of the synchronous detector 
from the overload. Second, the response time of each fre- 
quency band has been optimized. Theoretically, it is impos- 
sible to fine tune a filter as rapidly at 20 Hz as at 1 kHz, It 
simply takes a much longer time to detect a null at 20 Hz, 
but even more important is that the time response of a 20-Hz 
notch filter is much greater. In fact, a fundamental burst 



tone applied at the input of the notch will not immediately 
be nulled at the output. The entire input will first appear at 
the output and then decay away in proportion to f(notch)Q 
where Q is the Q of the notch, Thus, as the notch frequency 
is increased, so is the speed of the fine-tuning circuitry. The 
bottom frequency band (20-200 Hz) is relatively slow but 
nulls optimally at 20 Hz. For higher frequencies the instru- 
ment response time improves. 

Output Amplifier 

For the amplifier following the notch Filter, gain accuracy 
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Fig. 2. 8903 A notch filter fine tuning and balance circuits. 



fine-baJance circuit insures thai |V BP /V, N | - 1 at the input frequency 
(no amplitude error). The tine-tune circuit is described in some detail 
here; the fine-balance circuit operates in a similar manner. Fig. 2 is a 
simplified schematic of both circuits. 

The fine-tune circuit operates as follows (see Fig 2) The low-pass 
filter voltage V p, which is phase-shifted 90° from the input, drives the 
tune comparator USA, which turns JFET switch G5 on and off. When 
05 is on, point A is essentially grounded and no current flows into the 
integrator capacitor C3 When 05 is off, the notch amplifier feeds the 
tune integrator Because of this chopping action, the dc current that 
flows into the integrator is caused by the notch output component 
(hat is synchronous with V LP .* 

The tune integrator output is a dc voltage that can be changed only 
by the dc current flowing into the integrator. This voltage feeds one 
<nput of me multiplier U9A The other input to the multiplier is V fiP . The 
product of the two inputs is summed into the state- variable filter via 
U3 (see Fig. 1). The net result is that the tune integrator voltage can 
change the effective resistance (R T in Fig. t) and hence the center 
frequency of the notch filter 

The direction of change is such that any notch output components 
in synchronism with V lP are reduced. At steady state, the 

" Actually, th\i{i Item 'Olcft output can atso cause dc current to II-'.. .■ 

integrate* This error causes a snraH amounl of the fundamentaF frequency component to 
• r ?ijgn i he notch fitter. However tne maximum: error contribution re dsfortfon mea- 
&u?amerts ,mm 1hl£ offett is only about 0*6 <JS 



and good frequency response are important. In ac level 
measurements the signal travels through the output 
amplifier chain and is detected by the output rms detector. 
Any error degrades the H9Q3A's performance. The output 
amplifier allows the low- lev el dislorliun products leaving 
the notch to be accurately detected by limiting the required 
dynamic range of therms-to-ric converter. The rms detect nr 
is accurate over only a 30-dB range, and the output 
amplifier boosts these signals into the range of the detector. 
As with the input circuitry, rapid recovery from overload 
conditions is crucial. If the notch becomes rnistuned be- 
cause of a disturbance at the input, the output suddenly 



tune voltage is stable, although not necessarily zero, and no dc 
current flows into the tune integrator. Therefore, there is no notch 
output component synchronous with 

The fine- balance crrcuit insures that there are no components 
synchronous with V BP in the notch output V 0P and V LP are both at the 
fundamental input frequency and are in quadrature with each other. 
Thus, at steady state, there is no fundamental frequency in the notch 
output. 
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increases dramatically, sending the amplifier into overload. 
This in turn generates dc offsets in the amplifier chain that 
can take seconds to decay, even if the amplifier may start 
operating again much sooner. Thus a large low-frequency 
impulse appears at the amplifier output along with the 
signal being amplified. The ite signal is transferred 

to the output rms detector, which responds to a time- 
weighted average &i the total, lithe impulse is a significant 
fraction of the signal, the rms detector will not give a true 
indication of the signal amplitude. This can cause problems 
when the instrument ranges automatically- In effect, it 
forces the leveling algorithm to wait much longer In 1:011- 



AUGUST 1980 HEWLETT-PACKARD JOURNAL 15 



)Copr. 1949-1998 Hewlett-Packard Co. 



firm that the amplifier output is within the detector's range. 
This in turn slows down the rate at which the instrument 
can determine the proper measurement range and display a 
reading. Even if the impulse is 20 dB less than the signal, 
the detector error can be as much as 0,5%, This will also 
increase the amount of time required by the instrument to 
make an accurate reading from the detector once the proper 
range has been obtained. 

To alleviate these problems, the size of the transient is 
minimized in three ways. First, operational amplifiers and 
circuit configurations are used that have better than average 
overload immunity. Second, an overload detector is placed 
at the output of the rms detector, If the signal level becomes 
too large, the overload detector trips and the amplifier gain 
is reset to dB, Third, a 23-Hz high-pass filter is placed 
before the output detector. This significantly reduces the 
duration and amplitude of any transient and hence keeps 
the transient from significantly increasing the total mea- 
surement time. The only delay factors that remain are the 
controllable and predictable settling times of the notch 
circuit and the rms detector 

The response time of the output rms detector is a com- 
promise between rapid settling and low-frequency accu- 
racy. A configuration was selected that settles to within 1% 
of a 10:1 step in 350 milliseconds, and has a steady-state 
errnrofU.2%at 20 Hz, This includes filtering in the detector 
and additional filtering following the detector to reduce 
excess ripple. For leveling purposes the ripple is not sig- 
nificants so the microprocessor uses the detector out put 
when leveling and avoids the extra delay contributed by the 
additional filter. The output detector and filters could have 
been designed with switchable time constants to respond 
more rapidly for higher- frequency signals. However, the 
penalties would have heen additional circuit complexity 
and the ambiguity of not knowing when to invoke the longer 
time constant. A 20-kHz signal, for example, might still 
have a significant low-frequency component, which would 
cause excessive error with a more rapid time constant. 

Oscillator 

\I,;riY tA the design considerations for the notch filter 
apply equally to the oscillator. In both cases tuning consid- 
erations were the same, with switch able resistor networks 
used as the decade tuning elements and four-ijuadrant 
multipliers used for amplitude control. In many ways the 
oscillator and notch circuits can be seen as duals. The notch 
generates a pair of zeros on the jcu axis that reject the funda- 
mental component, while the oscillator generates a pair of 
poles on the j^axis that generate sustained oscillations. The 
trick in the oscillator is to keep the poles exactly on the axis 
to maintain constant output amplitude. This must be done 
continuously by the automatic leveling circuit. If the fre- 
quency of the circuit deviates from the desired frequency, 
the circuit can be fine tuned by the microprocessor, which 
monitors the output frequency on a sampled basis. The 
major performance goals of the oscillator were low noise 
and distortion, rapid amplitude and frequency settling, and 
digitally programmable frequency control- Again the 
state- variable filter configuration a long with a special level- 
ing circuit offered the flexibility and performance required. 
The oscillator design is described on page 10, 



It was determined during development that the oscillator 
would have to run at a constant output level to maintain 
reasonable settling and noise performance. It was also de- 
sired to have a floating output. The attenuator and output 
amplifier circuit (see page 12) lakes the oscillator output 
level and translates it to the selected floating output 
amplitude. To minimize cost and achieve overall output 
accuracy goals the attenuation is done in two stages. Coarse 
amplitude steps are implemented with a 2,5-dB. |l step at- 
tenuator network, Smaller steps are provided by a resistive 
ladder network that adjusts the amplitude linearly in small 
discrete steps, The combination can adjust the amplitude 
within a nominal ±0, 15% worst case. Computation of the 
proper switch settings is an easy job for the computational 
skills of the microprocessor. 

EMI Design 

Meeting the required electromagnetic interference (EMI) 
and susceptibility goals was a bit more challenging than 
initially expected. Large-amplitude RF fields tend to gener- 
ate voltages on exposed cabling and circuits. These voltages 
overdrive many of the active circuits, causing nonlinear 
operation and distortion. To avoid direct exposure to these 
fields, the analog circuits are housed in an internal EMI- 
tightbox. The box has an aluminum frame around the sides. 
The bottom cover is the ground plane of a printed circuit 
board and the top cover is a removable EMI- tight lid. Re- 
moval of the lid. which is held in place by only two screws, 
makes all the circuits available for service. The micro- 
processor boards are sufficiently shielded by the instru- 
ment cabinet and do not require the extra shielding. To keep 
the HP fields from developing voltages on the cabling feed- 
ing the circuits, special precautions were taken. First, from 
the inner box to the front panel, shielded cable is used. 
Second, BMC connectors are provided on the front panel. 
The BNC connectors allow the attachment of shielded ca- 
bles directly to the instrument if desired, thus preventing 
EMI pickup. The instrument's digital circuitry also gener- 
ates EMI related to harmonics of its 2-MHz clock, This 
problem was minimized by means of RF gaskets on some of 
the cabinet seams and by installing an EMI suppressing 
filter on the power line input. As a result, the instrument 
will not disturb sensitive receivers operating nearby, and 
yet will perform well near a powerful transmitter. 

Frequency Measurement 

A key feature of the 8903A is its ability to measure fre- 
quency automatically, even when the input waveform may 
have a significant amount of noise and distortion and the 
amplitude may vary from 6 mV to 300V. Part of this problem 
is solved because the instrument is autorangtng and keeps 
the leveled waveform within fi dB over most of the input 
amplitude range. But before the signal can be accurately 
counted it must first be converted into a binary signal hav- 
iriLs the same period as the major frequency component in 
the waveform, and herein lies a problem. If a zero-crossing 
circuit is used, noise may cause multiple crossings and a 
false indication of the frequency. Hysteresis in the detector 
will help, but if the hysteresis is too large, smaller- 
amplitude waveforms may not trigger the detector at all 
while large-amplitude waveforms will have relatively little 
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hysteresis protection when large noise components are 
present. To alleviate this problem* the 8903 A employs vari- 
able hysteresis. As the peak amplitude of the signal varies, 
so does the hysteresis level, which is maintained at approx- 
imately one-half the positive peak for the positive portion of 
the waveform and approximately one-half the negative 
peak for the negative portion of the waveform. Hence noise 
immunity remains constant regardless of the incoming 
waveform, Hysteresis is implemented with a bipolar peak 
detector followed by a dual comparator. The waveform is 
transferred to a reciprocal counter, which measures the 
period of the signal, and the microprocessor inverts this 
period to get frequency. 
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A Custom LSI Approach to a Personal 
Computer 



by Todd R, Lynch 



THE NEW IIF-85 PERSONAL COMPUTER, which 
was featured in last month's issue, is a system totally 
integrated into a single package. Included in this 
system are a CRT, printer, tape drive, and keyboard- To 
control the L'O (in put/output J devices and to interface to 
KOM (read-only memory) and RAM (random-access mem- 
ory), the design uses LSI (large-scale integrated) circuits 
designed and fabricated by Hewlett-Packard. 

The design of custom chips minimizes the cost of the 
electronics. Also, the power dissipation is reduced to a level 
ih^t permits the use of air-convection cooling, eliminating 
tlie need for a fan. The LSI designs save large amounts of 
printed-circuit-board space, making a small system pack- 
age possible. By designing LSI circuits dedicated to each 
mechanical subassembly, features can be added to the over- 
all system that would be nearly impossible to incorporate 
with discrete logic designs. 

HP-85 LSI Chip Set 

There are nine custom LSI circuits in the HP-85. These 
circuits are interconnected as shown in Fig. 1. Eight of the 
LSI circuit designs use NMOS technology, The ninth circuit 
uses a bipolar technology with two layers of metallization. 

The heavy emphasis on LSI design makes the HP-85 very 
compact. The electronics in the machine consists of these 
circuits plus a power supply, clock generator, CRT raster 
scan circuitry, and various motor and print head drive cir- 
cuits. 

The system Is partitioned as follows. 
* CPU (central processing unit]. The CPU commands the 

bus control lines and interfaces direct Iv to the system 



operating commands (firmware) stored in the four ROMs. 
ROMs, Kach ROM chip has the same basic custom design, 
differing only in the bit pattern permanently stored on 
each chip. 

Read,' Write Memory, The user memory consists of eight 
commercially-available dynamic 16KX1 RAMs. 
RAM Controller. This chip is designed to interface be- 
tween the CPU and the RAM*. 

Buffer. This chip is designed to provide the capability for 
system expansion by adding plug-in units in the rear of 
the machine. 
Keyboard Controller. 
Printer Controller. 
CRT Controller. 

Display Memory. Four commercial dynamic 16Kxi 
RAMs are used to store data for the display- Approxi- 
mately one-fourth of this memory is used for the alpha- 
numeric data and the remainder is used for graphic data, 
■ Cartridge Controller, 

Sense Amplifier. This bipolar chip interfaces the tape 
cartridge controller to the magnetic tape head. 
Table I is a summary of the custom LSI chip sizes, number 
of pins {wiring connection points), and power dissipation. 
The machines sixteen-hit address allows direct access to 
65.536 bytes of information. The memory map in Fig. 2 
shows how this space is allocated. Note that 32K bytes of the 
address space are devoted to ROM while 16K bytes are for 
RAM. Additional RAM and ROM capacity may he added to 
the system. The L'O address space occupies the upper 256 
bytes of memory. Each I/O device controller has from two to 
four dedicated addresses assigned to it, 
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Fig. 1 . The HP-85 system block 
diagram contains nine custom LSI 
circuit designs. The user and dis- 
play memories are the on fy parts tn 
the diagram that use commercially 
available circuits. The develop- 
ment of custom circuits enabled 
the system fo be contained in a 
compact package at low cost and 
without the need for a cooling fan. 
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Table I 






Custom LSI circuit characteristics All 


oi the 


s are fabricated 


with N'MOS silicon gale technology except for tlie sense amplifier 


which is made with dual-layer-metallization bipolar technology. 


Circuit: 


Size (mm) 


#Pins 


Power (mW) 




4 m 




330 


ROM 


5.41 


2fi 


200 


RAM Controller 


356 


1 


220 


Buffer 


2.54X3.35 


40 


340 


Keyboard Controller 


3,78x4.39 


42 


200 


Printer Controller 


4.7HX5.44 


40 


300 


CRT Controller 


4.14x5,51 


4C 


200 


Cartridge Controller 


3.63x3,86 


28 


55 


ose Amplifier 


1.50x1,55 


16 


150 



CPU Design 

Several commercially available CPUs were considered at 
the beginning of the project, but none could provide all the 
features needed to efficiently implement a powerful scien- 
tific BASIC language machine, BASIC requires high- 
precision arithmetic, which Is best accomplished with dec- 
imal rather than binary numbers, Mirny different stacks are 
needed to parse (separate a statement into executable steps) 
and execute the language. The ability to handle variable- 
length data is required for variable-length tokens (bit se- 
quences from one to several bytes in length), and multilevel 
vectored interrupts are needed lo handle the HO devices in 
real time. The design of the custom NMOS CPU incorpo- 
rates all of these requirements plus many more. 

Fig. 3 shows a block diagram of ibi-- HP-85 CPU. Major 
blocks are the 7000-bit PLA [programmable logic array], 
64-byte register bank, and eight-bit ALU (arithmetic logic 
unit) and shifter. Each CPU instruction is decoded by a 
microprogram in the PLA that directs the rest of the chip to 
perform the desired function, 

A very powerful feature of the CPU. which is incorpo- 
rated into the PLA microprogram, is [he ability to handle 
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Fig. 2. The HP8S memory map allocates fourBKbyte system 
ROMs and one 16Kbyte interna! RAM The RAM is used for 
user memory and 256 dedicated tiO locations. By using bank 
selection, up to six more 8K-byte ROMs may be added The 
RAM may be expanded by 16K bytes it a piug-m module is 
added to the system. 




Fig. 3. The CPU m the HP-85 is a custom NMOS circuit de- 
sign A 7K-pit programmable logic array (PLA) controls the 
interactions between the register bank, the ALU, and the rest 
of trie HP-85 system. 

data ranging from one to eight bytes in length. This multi- 
byte feature lets the programmer, for example, add two 
eight-byte mantissas, increment a sixteen-btt address, or 
load into the CPU a three- byte token from memory , each 
with a single instruction. With an off-the-shelf CPU this 
would require setup and iteration within a software loop. 

The eight-bit ALU is capable of shifting and can do both 
decimal and binary arithmetic. The program rner sets a 
mode bit to specify the type of shifting or arithmetic de- 
sired. This, combined with the multibyte feature, lets the 
programmer easily work with signed, floating-point man- 
tissas up to sixteen digits in length or two's- complement 
binary integers up In sixty-tour bits In length. 

Of the f>4 eight-bit registers contained in the CPU, one 
pair is dedicated to the program counter, one pair to the 
stack pointer, and one pair to internal index calculations. 
The rest are general-purpose registers. One advantage erf 
having a register as the program counter is that it can easily 
serve as one of the operands for any CPU instruction, If it is 
modified by an instruction, then an immediate jump to the 
new location occurs. The stack pointer points to the sub- 
routine return address stack stored in memory. Additional 
stacks can be created with any other consecutive pair of 
registers in the CPU. Thus, the programmer can maintain 
many different stacks at any one time. Sixteen instructions 
are dedicated to manipulating data on the currently desig- 
nated stack. Any pair of consecutive registers can also he 
used as a sixteen-bit index register. This lets the program- 
mer Index into several data arrays at the same time. 

To handle real-time I/O devices, the CPU has multilevel 
vectored interrupt handling ability. Up to 127 interrupt 
vectors can be accommodated by the architecture. The CPU 
can also be halted by an external device. This lets that 
device control the system bus, so it can have direct memory 
,n c I.*.-, rl f high speeds. 

No accumulator is present in the CPU. This is made 
possible by the design of the 64-byte register bank. The 
registers are constructed as a two-read/one- write memory, 
This means that, at a given time, any two bytes can be read 
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Fig. 4. Independent control of each word line is used to 
achieve a cell that can he readout onto either bit tine. The RAM 
made up with these celts can be read two words at a time. 



from the memory and operated on in the ALU. The result is 
relumed to one of the accessed byte locations. This se- 
quence takes one processor cycle (1,6 /is). 

To design a two-read/one- write memory, a standard static 
RAM cell with a special variation is used (see Fig. 4). In a 
normal static RAM cell, a common word line enables both 
transmission gates between the bit lines and the cell, in a 
two-read cell, there must be two different word lines. Word 
line A enables complement information from one cell to be 
read out onto the Q-bit line while word line B enables true 
information from another or the same cell to be read out on 
the Q-bit line. 

The circuit design of such a cell must be done carefully. 
Since the bit lines are precharged before reading, it is pos- 
sible that by enabling a single word line, a cell could be 
made to flip (change stale) rather than be read. Inadvertent 
flipping is prevented by using a low voltage for the word 
line and a proper size ratio between the cell's pull-down 
transistor and the transmission gate. 

The static RAM cell's pull-up device also received close 
scrutiny during the design phase. A cell with a low-power 
depletion-load pull-up requires an area larger than 6400 
square micrometres (10 square mils) and a quiescent power 
of 12U microwatts. By using a polysilicon pull-up resistor, 
the RAM cell area is reduced by 40% and the power is 
reduced by 80%. The savings in area amounts to a reduction 
in size of more than 0.25 mm on each side of the die. The 
reduction in total power consumption is fifty milliwatts. 

To obtain these reductions, a polysilicon resistor process 
was developed for the CPU. An additional masking step is 
required to define the regions on the chip where the 
polysilicon layer is lightly doped, hence creating the poly- 
silicon resistors, The doping level was chosen for a sheet 
resistance of 10 7 ft/a An un doped layer could not be used 
because the sheet resistance was so high that the junction 
leakage currents at elevated temperatures resulted in un- 
welcome voltage drops across some resistors. It was also 
discovered that if the contact mask overlapped the polysili- 
con resistor area, the aluminum metallization could spread 
into the lightly doped resistor and lower its resistance. 



System Control and Timing 

Eight bus and three control lines leave the CPU. The 
eight-bit bus is used to time-multiplex a sixteervbit address, 
instructions, and multibyte data quantities, The three con- 
trol lines — LMA (load memory address), RD [read) and WK 
(write) — indicate to the system circuits what type of infor- 
mation is on the bus. 

When LMA is low, all chips in the system know that one 
byte of a two-byte address will be placed on the bus. The 
LMA signals always come in pairs since all addresses are 
sixteen bits, Each chip in the system reads this address and 
compares it to the address range to which the chip is sup- 
posed to respond. Setting RD low indicates that the CPU 
wants to read the contents of the address most recently sent 
out. When WR goes low, the CPU wants to write into the 
location most recently addressed. 

The multibyte feature is accomplished by sending out an 
address with two LMA signals followed by one to eight RDor 
WR signals. The circuit being read from, or written to, is 
expected to increment its memory address register every 
time it sees a RD or WR. The CPU sometimes can fetch 
consecutive instructions from memory by simply sending 
out additional RD commands, This speeds up the instruc- 
tion rate of the machine since sending out an address for 
every instruction and piece of data is not necessary. 

Four nonoverlapping phases (Fig, 5) with 200 ns width 
and 200 ns spacing clock the circuits in the HP-85. In one 
cycle the system can access a pread dressed memory loca- 
tion, read it Into the CPU, add it to a CPU register and store 
the result in the CPU register. An eight- byte add requires 
eight cycles plus the time to fetch the command, which is 
usually three cycles. 

Before any chip writes on the bos, the bus gets precharged 
to a logic one, Therefore, the circuit desiring to put a one on 
the bus must merely continue to hold the bus high. If it 
wants a zero, it must discharge the bus, In NMOS t echo oh 
ogy It is easier to discharge a bus than to bring it high. 
Consequently, the circuits were designed with large de- 
vices to discharge the bus and relatively small devices to 
maintain a high level on the bus. This minimizes the chip 
area required for each circuit's driver logic. 

The system bus has fewer spurious transitions using the 
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Fig, 5, HP-85 system timing and control The control lines are 
valid by the time the $12 clock pulse occurs. If a system chip is 
to be read, it must respond during the <£2 clock pulse. Informa- 
tion then enters the CPU, goes through the ALU, and is stored 
durtng the next tfeT clock pulse. When the CPU sends out 
addresses or data, they are valid by the time the &1 clock 
putse occurs. 



20 HEWLETT-PACKARD JOURNAL AUGUST 19BQ 



)Copr. 1949-1998 Hewlett-Packard Co. 



precharged scheme. Timing is such that the driving circuit 
is not enabled until its data is valid and ready to go on the 
bus. It is not possible to discharge the bus inadvertently 
when the chip wants a logic one to be output, 

ROM 

The NMOS system ROM is organized as an 8K by eight- 
bit array. Because this circuit must respond to multibvte 
transfers* its memory address register is designed to incre- 
ment as RD commands are given. The circuit can be enabled 
or disabled by a bank select command. When plug-in ROMs 
are used, it becomes necessary to selectively enable or dis- 
able ROMs in the address space from 24K to 32K (Fig. 2). 
Each ROM in this address space has a unique eight-bit bank 
number, and only one bank can be active at any given time. 

The core of the ROM consists of 64- Input, minimum 
geometry, NAND gate arrays. The principle of operation is 
to precharge the NAND stack from both ends. This requires 
600 ns. At the beginning of phase #12., one 64-device stack 
per bit of output is selected. If the stack discharges, the bus 
bit will be a zero, otherwise it will be a one. The transistors 



in the stack axe preprogrammed during fabrication with the 
depletion mask, A depletion transistor is a logic zero, an 
enhancement transistor is a logic one. A high voltage on 
either transistor causes conduction, while a low voltage 
causes conduction only through the depletion device. 

When trying to discharge a 64-device stack. 63 inputs are 
high and the selected bit is low. If a depletion device is at 
that location, the stack discharges, producing a zero. 

RAM Controller 

The NMOS RAM controller chip interfaces the CPU to 
eight 1 6K dynamic RAMs. Because it is a memory controller 
chip it has an incrementing memory address register. Also, 
the address space to which it responds is mask program- 
mable. Therefore, a RAM controller chip different from the 
one in the mainframe of the HP-85 is used for the 16K RAM 
plug-in. This chip knows not to respond to dedicated I/O 
addresses in the upper 256 bytes of memory. 

An important function of the RAM controller chip is 
refreshing the dynamic memory. An internal timer on the 
chip tells it when the next sequential location is due to be 



The HP-85 Software Development System 



by Nelson 

The HP-85 is based upon a new eight-bit custom processor. One of 
the first steps in the development of the machine was to provide a set 
of software development tools. The HP-65 software development 
system consists of an assembler, a hardware interface for the HP-85 
simulator, and software debug system (see Fig. 1). 

The assembler for the HP-B5 processor was designed to run on an 
HP 1000 Computer The assembler supports the full range of the 
HP-85 processor's instruction set and provides several useful 
pseudo-operations Programs may be either absolute or relocatable 
and the program origin may be reset at any time. Several instructions 
are provided to facilitate data definition and may exist locally withm 
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Fig. 1 . HP-85 software development system 



A. Mills 

the program, or may be retrieved from a file of global data definitions. 

The hardware interface designed for the HP-85 software develop- 
ment system provides the ability to use an HP 1000 to control execu- 
tion of the HP-85 simulator The interface contains an on-board RAM 
which can be downloaded from the HP 1000 and is used to simulate 
the HP-85 ROM Included are two breakpoint registers which can be 
used to halt execution of the HP-85 processor at either of two 
specified addresses. The processor can be made to execute steps in 
either a continuous-run mode, arm a single-step mode where execu- 
tion halts after each instruction. The interface also provides the ability 
to halt the HP-85 processor at any time, and to read data from or write 
data to the system RAM or CPU registers 

The software debug system, provided as part of the development 
system, runs in the HP 1000 and performs lour important functions. It 
includes a relocatable loader that is used to download (he RAM on the 
interface board with the software under development. It controls 
execution of the HP-85 simulator via the interface board, ft displays 
the current status of the breadboard, including CPU status, program 
counter, and the current contents of all CPU registers and any 
specked memory locations Finally, it provides the means for system 
developers to modify the contents of memory, registers, status, or 
program counter, and to clear and set breakpoints. 
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refreshed. The chip then waits for a cycle when memory is 
not being accessed to do the refresh. There is a two-element 
queue in case a vacant memory cycle does not come before 
the following location must be refreshed. 

Keyboard Controller 

The circuit needed just to interface to the keyboard would 
have been quite small, Consequently, four programmable 
timers were added to this NMOS chip. The timers can count 
up to 27 hours with one-millisecond resolution. Each has a 
maskable interrupt. When the timer gets to its preset count, 
it interrupts, resets and counts again. 

The keyboard portion of the chip can scan an 8xl0-key 
keyboard plus three dedicated keys — shift, control, and 
caps lock, The key debounce time is mask programmable 
from 1.67 to !Tti9 ms. After the key has been denounced, 
the chip generates an interrupt. An internal ROM is used to 
convert the key position to its ASCII* equivalent. 

Another feature of the chip is an output to a speaker for 
audio tones. Via firmware, a 1,2-kHz output tone can be 
obtained, or the frequency can be varied by periodically 
setting and clearing an internal flip-flop. This output can be 
accessed by the programmer using a command that can 
specify the frequency and duration of the tone- 
CRT Controller 

The NMOS CRT controller 1 interfaces the CPU to a 
127-mm diagonal CRT and its dedicated display memory. 
The memory is four 16K x 1 dynamic RAMs. This is enough 
storage to hold a full display of graphics in formation 
(25HX192 dots) plus four displays of alphanumeric infor- 
mation (4x32 characters 'line x 16 lines/display). The chip 
must make sure that this memory is properly refreshed 
during the vertical retrace. The interface to the video drive 
circuitry consists of a vertical sync pulse (60 Hz) , a horizon- 
tal sync pulse (15.7 kHz), and a video line (4,9 MHz). An 
internal ROM provides a dot pattern translation of the 
stored ASCIT characters. 

Printer Controller 

The NMOS printer controller 2 is the interface between 
the CPU and the thermal moving- head printer, To perform 
this function, the circuit must control a printhead drive 
motor, a paper advance motor and an eight- dot printhead. 
An internal 32-byte RAM allows firmware to buffer one 
print line of alphanumeric data at a time. An internal 
ROM translates an ASCII character into its appropriate 
dot pattern. 

Tape Controller and Read/Write Amplifier 

The NMOS tape controller and the bipolar read/write 
amplifier 1 control the tape unit in the HP-85. The read/write 
amplifier reads and writes data on the tape through a two- 
track magnetic head. The tape is formatted using a 1:1,75 
delta distance code, in which 8- kHz flux reversals represent 
zeros and 4.6-kHz flux reversals represent ones. 

The tape controller encodes digital information into this 
format and sends it to the read/write amplifier. When read- 
ing the tape, the signals from the read/write amplifier are 
decoded by the controller. The tape drive motor direction 
and speed are also controlled by the chip, 

"American Stendd'd Code lor Information rn tertian 34= 



Buffer 

The NMOS buffer is the ninth chip in the system. With its 
help, the system bus can be expanded to include all of the 
plug-in I/O slots. It is capable of driving a 150-pF load. The 
delay between NMOS- level input and NMOS-level output 
is 50 us. The signals that the buffer passes are eight bus 
lines, three control lines, and the interrupt and halt lines. 
The bus and control lines are designed to be bidirectional. 

The clock lines are not buffered since this would cause 
unwanted skewing. The clock generator is capable of driv- 
ing all internal as well as all external loads* 

Acknowledgments 

A great deal of credit must be given to the team that made 
this LSI chip set possible. Tim Williams designed the RAM 
controller and buffer. Donn VVahl designed the ROM. The 
keyboard controller was an effort from Jim Hutch ins and 
Jerry Erickson. ferry also designed the CRT controller. The 
printer controller was the work of Clement Lo. Implementa- 
tion of the NMOS CPU was done by Jim Axtell. The car- 
tridge controller design is that of Doug Collins. The sense 
amplifier is the result of the combined efforts of Mike 
Moore, Doug Collins and Mike Barbour. Special thanks 
must be given to Don Hale and Howard Shishido, who did 
the majority of the chip mask designs, Howard Honig con- 
tributed to the project by writing most of the test programs 
for the LSI circuits. Mike Pan. Rock Davidson, and Tom 
Kraemer also contributed to the chip set design effort. Bob 
Tillman and Steve Larsen helped get the NMOS process and 
polysilicon resistor development up to production stan- 
dards, Payne Freret contributed some good ideas to the CPU 
design. Rick Bell did the design of the clock oscillator and 
generator. 

Reference 

1, J.F. Bausch. 'A High-Quality CRT Display for a Portable Com- 
puter," Hewlett-Packard Journal, July I960. 

2, C,C Lo and R,W, Ketl, "A Compact Thermal Printer Designed 
for Integration into a Personal Computer. " Hewlett-Packard Jour- 
nal, July 1980. 

3, D.J. Collins and B.C. Spreadbury, M A Compact Tape Transport 
Subassembly Designed for Reliability and Low Cost," Hewlett- 
Packard Journal, July 19S0. 



Todd R. Lynch 

Todd Lynch is a native of Rochester, 
New York. He joined HP Laboratories in 
1 972 after working for a year on process 
control logic designs. Since coming to 
HP he has worked on computer ar- 
chitectures, and since transferring so 
what is now HP's Corvallis Division, 
Todd designed the architecture for the 
HP-85 CPU He was the project man- 
ager for the integrated circuits in the 
HP-85 and is now project manager for 
high-end mainframes. He received a 
BS degree from Grove Cay College, 
j* Pennsylvania in 1970 and an MSEE 

_^ft (rom the University of Wisconsin in 

1971 Todd is co-chairman for the Willamette Valley Junior Tennis 
Tournament and in addition to tennis enjoys woodworkjng, soft ball, 
and skiing. He lives in Albany, Oregon with his wife and new son 




22 HE'ALETT-PACKARP JOURNAL AUGUST 1980 



)Copr. 1949-1998 Hewlett-Packard Co. 



Handheld Calculator Evaluates Integrals 

The HP-34C is the first handheld calculator to have a key 
that performs numerical integration almost automatically, 
it may change your attitude towards what used to be 
regarded as a dreary tedious task. 

by William M. Kahan 



NUMERICAL INTEGRATION has been the subject 
of about two thousand books and learned papers, 
with a dozen or so "new" methods published every 
year. And yet the task in question has a simple geometrical 
interpretation seen in Fig. I: given an expression f(u] and 
lower and upper limits y and x respectively, the value 



I 



l*hu\du 



represents the area under the graph of f[uj for u between y 
and x. Why so much fuss? 

As I write this an electrical engineering colleague, Pro- 
fessor J. R, Woodyard, enters my office and asks to have 



Jo l u-1 



In u ' 



evaluated on my 1IF-34C Calculator (Fig, 2], Let's do it. 

Step 1. Key Into the calculator under, say, label A a pro- 
gram that accepts a value u in the display [X regis- 
ter) and displays instead the computed value of tin* 
integrand 



\ u(u 



l/(ln u) 



Step 2. 



Step i 



Step 4. 



Fig. 3 shows an IIP-34C program to do this. 
Restore ihe calculator to RUN mode and sei the dis 
play to, say, FIX 5 to display five decimal digits after 
the point, which are as many digits of the integrand 
as my client says he cares to see. fMore about Mi is 
later.) 

Key in the lower and upper limits of Integration 
thus, ENTER T 1, thereby putting into the Y reg- 
ister and 1 into X. 

Press jy A, wait 25 seconds until the display shows 
0,03662, then press x^y to see 0.00001. We 
have just calculated 



« What can go wrong and how do we avoid it? 

What else have we learned? 
These questions and others are addressed in the following 
pages, 

Tolerance and Uncertainty 

Inlegrals can almost never be calculaled precisely. How 
much error has to be tolerated? The Jy key answers this 
question in a surprisingly convenient way. Rather than be 
told how accurately 1 = /y f[u|du should be calculated, the 
HP-34C asks to be told how T many figures of fju| matter. In 
effect, the user is asked to specify the width of a ribbon 
drawn around the graph of flu), and to accept in place of 
I an estimate of the area under some unspecified graph lying 
entirely within thai ribbon. Of course, this estimate could 
vary by as much as the area of the ribbon, so the calculator 
estimates the area of the ribbon too. Then the user may 
conclude from Fig, 4 that 

/ = fnrea under a graph drawn in the ribbon) 
area of the ribbon \ 

The calculator puts the first area estimate in its X register 
and the second, the uncertainly, in the Y register. 

For example, ffu] might represent a physical v\\m.\ whose 
magnitude can be determined only to wilhin. say, i 0.005. 
Then the value calculated as f(u) is really f|u) ± Af(u) with 
an uncertainty iffu) = GJHJfj. Consequently FIX 2, which 
tells the calculator to display no more Hum two decimal 
digits alter the point, is used to tell the calculator that 
decimal digits beyond the second cannot matter. Therefore 
the calculator need not waste time estimating I ± 41 = 
Jy |f(u} ± Af(u))du more accurately than to within ail 
uncertainly M - |/*Af(u)du|. This uncertainty is 
mated together with J ± AI, thereby giving Ihe calculator's 
user ;i loir idea of the range of values within which / 
must lie, 



J, 



0JJ3BB2 ± U.nnUUl, 



That w r as easy — too easy. Woody ard smiles as if he know 
something 1 don't know. Could the calculator be wrong? 
How do es t h e ca I cu 1 a to r kn o\v the error 1 ies wi th i n ± 0. 00 1 ? 

Many other questions come to mind: 

Why is numerical integration impossible in general? 
i Why do we persist in trying to do it anyway? 

How do we do it? How well do we do it? 
i How does the Jy key compare with other integration 

schemes? 




u = y 



Fig. 1. An Integral interpreted as an area. 
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Fig. 2. HP-34C Calculator has keys to solve any equation and 
to compute integrals. 

The uncertainty Affu) is specified by the user via the dis- 
play setting. For instance. SCI 5 displays six significant 
decimal digits, implying that the seventh doesn't matter. 
The HP-34C allows the user's f- program to change the dis- 
play setting, thereby providing for uncertainties M(u) 
that vary with u in diverse ways, But users usually leave 
the display set to SCI 4 or FIX 4 without much further 
thought. 

By asking the user to specify Af[u] instead of Al the 
HP-34C helps avoid a common mistake — wishful think- 
ing. Other integration procedures, which conventionally 
expect the user to specify how r tiny A! should be, blithely 
produce estimates of J purporting to be as accurate as the 
user wishes even when the error Af(u] is far too big to justify 
such claims to accuracy. The HP-34C does not prevent us 
from declaring that f(u) is far more accurate that it really is, 
but our attention is directed to the right question and not dis- 
tracted by questions we cannot answer. Whether we spec- 
ify Af after a careful error analysis or just offer a guess, we 
get estimates I ± A/ that we can interpret more intelligently 
than if w r e got only I with no idea of its accuracy or inaccuracy. 

A Survey of Integration Schemes 

Students are taught the fundamental theorem of calculus: 



as a process, applied to expressions, that starts with f and 
ends with F. But in professional practice that process 
hardly ever succeeds. A compact expression F{u) is al- 
most always difficult or impossible to construct from any 
given f(u]> For instance, neither 

/_ exp( -Lr'2jdu VZtt nor J t>xp(-u +■ x In u]du 

possesses a closed form, that is, an expression involving I <i i J \ 
finitely many elementary operations [+, - , X, -r T In, exp, tan. 
arctan. ,..) upon the variable x. Nevertheless, both integrals 
can be approximated arbitrarily accurately by aptly chosen 
formulas. So often do statisticians and engineers need 
values of those integrals that formulas for them, accurate 
to ten significant decimal digits, can now be calculated in 
a few seconds by pressing a key on certain handheld cal- 
culators, (Press Q on the HP-32E to get the first integral, 
the cumulative normal distribution; press xl on the HP-34C 
to get the second integral, the gamma function Ffl+x), 
whether x be an integer or not,) 

Almost every rare integrand f(u) whose indefinite in- 
tegral F(x] = J * f(u)du •¥ c is expressible in a compact or 
closed form can be recognized by a computer program that 
accepts the string of characters that defines f and spews out 
another string that represents F, [Such a program is part of 
I he MACSYMA system, developed at MIT, that runs on a 
few large computers — tw T o million bytes of memory — at 
several universities and research labs,) Perhaps the terms 
" compact* ' and ''closed form" should not be attached to the 
expression F(x), since usually, except for problems as- 
signed to students by considerate teachers, the integral 

Begin with u in the X register 
Save u in the stack 



\ u 



Recall u 



...:-N U.'(U-l) 

Recall u again 
, . - 1;ln(u) 






)du = 



F(x)-F(yl provided ~ F(u) ■ 
du 



flu I 



This means that one could calculate f if one could dis- 
cover somehow an expression F(u) whose derivative is 
the given expression f(u). Students are taught integration 



Display \ u'(u-1) - i;in u 

Fig. 3. This program makes the HP-34C calculate the inte- 
grand Vu/fu-T) - Hin u when the argument u is in the X 
register and key A ts pressed. Labels S, 0, T , 2, or 3 would have 

served as well as A 
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Al - V* Area of Ribbon « J^* Af(u) du 



Area! = f f * f(u) du 



f[u)-Af(u) 



Fig. 4, fne graph of an uncertain integrand iiu - A . dr? run 
any where m the ribbon bounded by the dashed itnes . The area 
under such a graph, l±i/. is uncertain by -t&l t which is one- 
half the area of the ribbon The HP-34C displays Us estimate of 
I ± A/ tn its X regtster and holds an estimate of A/ m its Y register 

F(x) far exceeds the integrand f(u 1 in length and complex- 
ity. Shown in Fig; 5 are two compact forms and one closed 
form forF(x] when f(u)=l'( l+u^j. The extent to which K(xj 
is here more complicated than f(u) is atypical I y modest out 
of consideration for the typesetter. The formulas in Fig. 5 
will remind many readers of hours spent on calculus prob- 
lems hiil they do not provide economical ways to calculate 
F(x) for any but very big or very tiny values of x. When I use 
the HP-34C's Jy key to caJculate F(l) = Jo du (T+u 64 ) = 
0,989367 ± 0,000004 the answer appears in 2uU seconds 
including 20 seconds taken to enter the f-prngrarii plus 180 
seconds for a result (in SCI 5). Calculating F[l) from any 
formula in Fig. 5 takes at least about ten times longer, not 
including the time taken lo dedui e the formula. Kngineers 
and s< ientists have Ions been aware of the shortcomings oi 
integration in closed form and have turned to other 
methods. 

Perhaps the crudest \\<i\ to evaluate /£ F[u)du is to plot 
the graph ol I[lj|, like Fig, 1. on uniformly squared paper 
and then count the squares that lie inside the desired 
area. This method gives numerical Integration its other 
name: numerical quadrature Another way, suitable £01 
chemists, is to pint the graph on papei nf uniform density, 
1 ui hji the ciird in question, and weigh tt. Engineers used 
to measure plotted areas by means of Integrating engines 
called piani meters. These range from Inexpensive hatchet 
planimeters of low accurac} tn Mviss-niade museum 
pieces costing hundreds of dollars and capable of three 
significant decimals, (For more details see reference 1). 
Nowadays we reckon that the computer will drive 
the graph plotter so it might as well integrate I no, 

Today's numerical integration techniques are best ex- 
plained in terms of averages like 

A~ J/|x-y) -/ v v Ffu}duy(x-yj 

which is tailed "the uniformly weighted average of f|u) 
over the interval between xand y ." Another kind of average. 



A = 2*w ( f(U|J where w,>o and Mv, 1 
I i 

finite weighted average of n samples itujh f|u 2 ) l'[u n l 



Provided the sample arguments u v u 2 * .... u n , called nodes, 
all lie between x and y the sample average A will approxi- 
mate, perhaps poorly, the desired average A^ and hence 
provide I = (x — y)A as an approximation to i = (x-yjA. 
Statisticians might be tempted to sprinkle the nodes u^ ran- 
domly between x and y — that is what Monte Carlo methods 
do. But randomness is a poor substitute for skill because 
the error A- A tends to diminish like 1 \ n as the number 
n of random samples is increased, whereas uniformly 
spaced and weighted samples provide an error that dimin- 
ishes like 1 n~. Other more artful methods do even better. 
Different numerical integration methods differ princi- 
pally in the ways they choose their weights Wj and nodes Uj. 
but almost all have the following characteristics in com- 
mon. Each average A is associated with a partition of the 
range of integration into panels as shown in I ig. B. Each 
panel contains one node Uj whose respective weight is 

Wj = (width of panel j)'( width of range of integration). 

The formula given above for A amounts to approximating 
the area in each panel under the graph of f|u] by the area of 

I angle as wide as the panel and as high as the sample 
tin, 1. The simplest method is the midpoint rule, whose 
nodes all lie in the middles of panels all of the same width. 
Other methods, like the trapezoidal rule and Simpsons 
rule, vary the panel widths (weights) and nodes in ways 
designed to exploit various presumed properties of the 
integrand f(uj for higher ai Which method is best? 

If this question had a simple answer there would not be so 
many methods nor would we need texts like "Methods of 
Numerical Integration" by IV I . Davis and P. Rabi now it/,- 
which contains Hi FORTRAN programs and three bib- 
liographies with well over 101 HI citations. 

For example, consider Gaussian quadrature. This method 

idely regarded as "best"' in the sense that it very often 
requires fewer samples than most other methods to achieve 
an average A that approximates the desired A to within 
some predesigned tolerance. Hut the weights and nodes of 

sian quadrature take quite a while lo calculate. Pro- 
grams to do so. and the resulting tables ol w eights and nodes 
for various sample counts n. have been published.-* Had 
we chosen Gaussian quadrature ha the |y key we would 

F{xjsf f(u)du where f(u) VM+u 64 ). 



Fixi *2,( 



-x 6 V.M64k + 1) jf ** 1 



^ esc (^) stgn tx) + x V ( x *V t (64k 1) if x : * 1 



; | [sin* «dM'(%^.+ ttmfti In fl i — ^ ]\ 

\ 2cos# k 7 

* (^ csc (^|) sign (x) \l tf 1) where tt k - (k-V 2 )7n32 



Fig, 5. Format integration transforms many a simple expres- 
sion t{u) into messy formulas F(x) of limited numerical utility 
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Fig. 6. The integral, regarded as an area, <$ here divided into 
four panels each of whose areas is approximated by The area 
of a rectangfe as wide as the panel and as high as a sample. 

have had to store a! least as many nodes and weights as we 
could expect to need fof difficult integrals, amounting to 
at least several hundred 13-digil numbers, in read-only 
memory. Bui that would have left no space in the HP- A4C 
for anything else, so a different method had to be found 

The Jy key could not use a method that generates just 
one average A because that gives no indication of how 
accurately it approximates A. Instead we looked only at 
methods that sample repeatedly and with increasing 
sample counts n l < n 2 < n a < ... to produce a sequence of 

increasingly accurate averages A T . A z > Aj Provided 

that sequence converges to A so fast that each \A^ X - A \ 
is considerably smaller than its predecessor, the error 
|A k -^ can be approximated accurately enough by 
| A^-Afc n | . ami the last average A^-i can be accepted in 
lieu of A as soon as | A^-A^.^ | is tolerably small, 

How small is "tolerably small"/ That depends upon the 
area of the ribbon discussed above under "Tolerance and 
Uncertainty.*" Since the integral f - Jy f(u)du inherits an 
uncertainty AI = | J~£ Af(u)clu| from the uncertainty if(u] 
in the integrand, so does A - f.'jx-y) inherit an uncertainty 
AA = AJ/|x-y| T which may be approximated by 



AA 



1=1 



in the same w r ay as A is approximated by A. Indeed. A and 
AA can be computed together since they use identical 
weights and nodes. And so the sequence A 1( A 2 * A s . ... is 
accompanied by a sequence of respective uncertainty esti- 
mates AA^ AA 2 . \A :i Now "tolerably small" can be 

defined to mean "rather smaller than AA k + 1 /" 

The foregoing argument provides an excuse For accept- 
ing Afc+i in lieu of A whenever two consecutive estimates 
A k and A^ tl agree to within AA k+1 , but it provides no 
defense against the possibility that convergence is not so 
fast, in which case A& and Aj^j might agree by accident 
and yet be both quite different from A, The Jy key waits for 
three consecutive estimates A k , A k+ll and A^+2 t° agree 
within AA kf2 - Only the most conservative integration 
schemes wait that long. While this conservatism strongly 
attenuates the risk of accidental premature acceptance of 
an estimate, the risk that three consecutive estimates might 
agree within the tolerance and yet be quite wrong cannot 
be eliminated. Later, under 'How to Deceive Every Nu- 



merical Integration Procedure," some such risk will be 
proved unavoidable, but the risk now is so small that 
further attenuation is not worth its cost 

The combination of ignorance with conservatism is 
surprisingly costly. Had we known in advance that A^ 
would be accurate enough we would have calculated none 
of the other averages. Instead, waiting for three consecu- 
tive averages to agree could easily cost some methods al- 
most fj.25 times as many samples as if only A^ had to be 
calculated, and more than I hat if the sample counts n T > 
n 2t n 3 , ... are not chosen optimally. For the/y key we chose 
nfc = 2 -1 and we chose a method whose successive aver- 
ages each share almost half of the previous average's 
samples, thereby preveoting the cost of ignorance from 
much exceeding a factor of 4. 

Memory limitations precluded the use of another family 
of methods known as adaptive quadrature. These methods 
attempt to distribute nodes more densely where the inte- 
grand f(u) appears to fluctuate rapidly, less densely else- 
where where f[u) appears to be nearly constant or relatively 
negligible- They succeed often enough that the best 
general -purpose integrators on large computers are adap- 
tive programs like Carl de Boor's CADRE; this and others 
are described in reference 2. Alas, adaptive programs con- 
sume rather more memory for scratch space than the 
twenty registers available in the HP-34C 

What Method Underlies the/* Key? 

The IIP-34C uses a Romberg method: for details consult 
reference 2. Several refioements were found necessary. In- 
stead of uniformly spaced nodes, which can induce a kind 
of resonance or aliasing that produces misleading results 
when the integrand is periodic, the Jy key's nodes are 
spaced noiiunifcirmly* Their spaciog can be explained by 
substituting, say, 

3 



into 



1= J f(lljdll= f f(~L.- V) I (l-lT)dv 

and distributing nodes uniformly in the second integral. 
Besides suppressing resonance, the substitution confers 
tw r o more benefits, One is that no sample need be drawn 
from either end of the interval of integration, except when 
the interval is so narrow that no other possibilities are 
available, and consequently an integral like 



\* 



J" 

Jn 



SIM du 

u 



won't hang up on division by zero at an endpoint. Second. 
1 = jy i(ujdo can be calculated efficiently when 
f(u] = g(u]V |x-u| or g[u]V[x-u](u-y] where g(u) is 
everywhere a smooth function, without any of the expedients 
that would otherwise be required to cope with the infinite 
values taken by the derivative f'(u] at u = x or u = y. Such 
integrals are encountered often during calculations of areas 
enclosed by smooth closed curves, For example, the area of 
a circle of radius 1 is 
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Jg \ u(4-u)du = 3.14159 ± 8.8x10' 

which consumes only 60 seconds when evaluated in SCI 5 
and only 110 seconds to get 3.141592654 = 1 4 - 10" 9 
in SCI 9. 

Another refinement is the use of extended precision. 
13 significant decimal digits, to accumulate the sums that 
define A^, thereby allowing thousands of samples to be 
accumulated, if necessary, without losing to roundoff any 
more information than is lost within the user's own f- pro- 
gram. The last example's 10 significant decimal digits of 
*r could not have been achieved without such a refinement. 

How Does the J J Key Compare with Other Integrators? 

What most distinguishes the HP-34C T s jy key from all 
other schemes is its ease of use. No step-size parameters, no 
plethora of error tolerances, no warning indicators that 
"can usually be ignored." Only the minimum informa- 
tion needed to specify f J (f(u) *Af[u])du has to be supplied. 
And because the /£ key is effective over so wide a range of 
integrals it ranks among the most reliable procedures avail- 
able anywhere- Usually it is tar faster than simpler proce- 
dures like the trapezoidal rule or Simpson's rule com- 
monly used previously on calculators. For integrands de- 
fined by programs that fit comfortably into a mid-sized 
handheld calculator that can hold at most 210 program 
steps, the Jy key is comparable in speed (count the number 
oJ samples] with the integrators available on large com- 
puters. For much more complicated integrands the best 
adaptive integrators on large computers are appreciably 
faster. 

One of the HP*34C's most important components is its 



Owner's Handbook. It is for most owners the first guide to 
the foothills of an awesome range of new possibilities. Two 
chapters are devoted to Jy. The first is introductory, and 
allows the user to evaluate simple integrals effortlessly 
and confidently. The second chapter is a longer explana- 
tion of the power and the pitfalls, concerned mainly with 
numerical integration generally rather than with the 
HP-34C in particular, This chapter had to be included be- 
cause its explanations and practical advice are not yet to be 
found in any text likely to be consulted by an owner, nor 
are they supplied by the instructions that accompany other 
integrators on other computers or calculators. This second 
chapter is part of the educational burden that must be 
borne by innovators and pioneers. The Owner's Handbook 
provides no formulas for the nodes and weights used by 
the HP-34C because they are not needed to understand 
how the Jy key works; instead they can be deduced from 
information in this article. 

Every numerical integrator like jy, which executes a 
user-supplied program to get the integrand's value f{u), 
imposes constraints upon that program. Some constraints, 
like requiring f to have a smooth graph on the interval of 
integration, are practically unavoidable, Others are nui- 
sances like 

■ Begin the f* program with a special label, say A'. 

■ Do not use certain memory registers . say #0 - #5. 

■ Do not use certain operations, say = and CLR. 

The Jy key is encumbered with no such nuisances. Thr 
f-program may begin with any of several labels, so several 
different integrals can be calculated during one long com- 
putation. The t-program may use memory registers freely 
and may use any operation key except Jy itself. (Due oi 



Evaluation in RUN Mode 



integrand in PRGM Mode 



Equation in PRGM Mode 



. get u 



. tor v 



save "', get u 



.lifv 



Wait for answer / 
to be displayed 



1 +mi* m 



, i'-u + ln(l+ur>e u ) 



Fig. 7. A program to evaluate 
i = f udu v(u) where v v\ ■ 
satisfies v-tt+ln (1 tuve J )=G 
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those keys is the HP-MiTs powerful SOLVE key* 4 Conse- 
quently this calculator is currently the only one that can 
evaluate conveniently integrals of implicit functions, 
Far example, let v = v{u) he the root of the equation 



ln(l 



Then 



/J U du!v[u) = 1.81300 ± 0,000005 



results from a program rather shorter than on any previous 
calculator; it is exhibited in Fig. 7. 

Furthermore, Jy may be invoked, ]ike any other function, 
from within a program, thereby permitting the HP-34C to 
SOLVE equations involving integrals, For example, 
solving 

J* cos [x sin 6)dti = 

for x - 2,405... takes a short program contained in the 
Owner's Handbook, and exhibits the first zero of the Bessel 
function Jnfx). 

How to Deceive Every Numerical Integration Procedure 
Such a procedure must be a computer program — call it 
P — that accepts as data two numerical values x and y and 
a program that calculates f(uj for any given value u T and 
from that data P must estimate J = /* f(u] du. The integra- 
tion procedure P is not allowed to read and understand 
the f- prog ram but merely to execute it finitely often, 
as often as P likes, with any arguments u that P chooses. 
What follows is a scheme to deceive P. 

First ask P to estimate / for any two different values x and 
y and forf(u)^0. Record the distinct arguments U| ■ u 2 , ,,.. u n 
at which P evaluates f(u). Presumably when P finds that 
ffu,] = f[u 2 ) = ... = f(u n J = it will decide that 1 = Q and 
say so. Next give P a new task with the same limits x and y 
as before but with a different integrand 

f(ll] = [(u-lhWli-UjiJ-^-l'U-UnH-. 

Once again P will calculate f(ui). f(u^) and finding no 

difference between the new f and the old, P will repeat 
exactly what it did before. But the new integral i is quite 
different from the old, so P must be deceived. 

The HP~34G + s jy key can be hoodwinked that way, Try 
to evaluate/-*^ f(u)du using first f(u) = programmed 
in a way that pauses fuse the PSE key] to display its argu- 
ment u, The calculator will display each sample argument 
it uses, namely f). ±HB. ±47 and *117, Next program 

f(u) * [u(u-88)(u + 88)(u-47)(u-h47][u~117)(u + il7)) 2 

and evaluate J_J| a f(u]du again, The calculator will say 
that both integrals are 0, but the second polynomial" 1 -; 
integral is really 1,310269 X 10 . That polynomial's graph, 
shown in Fig. 8, has the sharp spikes that characterize 
integrands troublesome for every numerical integration 
procedure. To Calculate the integral correctly* reevaluate 
it as 2 jo 12 f(ujdu t thereby doubling the spikes' width com- 
pared with the range of integration. 



The threat of deceit impales the designer of a numerical 
integrator upon the horns of a dilemma. We all want our inte- 
grators to work fast, especially when the integrand !{uj 
is very smooth and simple like f(u) = 3u— 4. But if the inte- 
grator is too fast it must be easy to deceive; fast integration 
means few samples f(Uj). implying wide gaps between some 
samples, which leave room for deceitful misbehavior. 
Figs. 9a- 9e illustrate the dilemma with two estimates of 
Jy f(u]du. The first estimate is based upon the three sam- 
ples drawn at the white dots, the second upon seven sam- 
ples including those three white plus four more black dots, 
fig. 9a shows why all sufficiently smooth graphs f(u) that 
agree at all seven samples have nearly the same integrals, 
but Fig. 9b shows how two integrands could provide the 
same samples and yet very different integrals. The coinci- 
dence in Fig. 9b is unlikely; successive estimates based 
upon increasingly dense sampling normally would reveal 
the difference as in Fig. 9c. However, situations like those 
illustrated in Figs. 9d and 9e are very likely to deceive, 

Textbooks tell us how T to avoid being deceived: avoid 
integrands f(u) among whose first several derivatives are 
some that take wildly different values at different places 
in the range of integration. Or avoid integrands f[u) that 
take wildly different values when evaluated at complex 
arguments in some neighborhood of the range of integra- 
tion. And if wild integrands cannot be avoided they must 
be tamed. We shall rejoin this train of thought later, 

Improper and Nearly Improper Integrals 

An improper integral is one that involves w in at least 
one of the following ways: 

One or both limits of integration are ±^, e.g., 




117 128 



Fig. S, The polynomial f(u) was devised to deceive the 
HP-34C into miscalculating its integral as instead of 
7 3i x1Q? a . This sptky graph ts typical of integrands that can 
baffle any numerical integrator. 73% of the area under air- 
graph lies under two spikes whose widths span fess than 4 % of 
the area of integration 
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* Trul 



(aj 



[b) 



(c) 



,d) 




M 



Fig, 9, Few samples (open circles) mean fast integration but a 
large possibility of error. More samples (soiid dots plus open 
circles) usually mean more accuracy, but not always, as in (h), 
(d) t and (e). (a) Which is the graph of Hup No matter, both 
have almost the same integral, (b) Which is the graph of f(u)7 
They have very different integrals, (c) Here two graphs that 
coincide on the first samples Q are distinguished by a signifi- 
cantly different outcome after second samples • are drawn. \d) 
If the graph of f{u} has a few sharp and narrow spikes, they will 
probably be overlooked during the estimation ot the integral 
based on finitely many samples, (e) if the graph of f(u) has a 
Step that was not made known during the estimation of the 
integral, then the estimate may be mistaken 



J_ r exp (-u") du = \ -2, 



The integrand tends to ±^ 
integration. e,g,. 



someplace in the range of 



ii(u)du = L 

■ The integrand oscillates infinitely rapidly somewhere in 
the range of integration, e.g., /q cos (In u]du = 

Improper integrals are obviously troublesome. Equally 

troublesome, and therefore entitled to be called nearly 

improper, are integrals afflicted with the foil owing malady: 
The integrand or its first derivative changes wildly with- 
in a relatively narrow subinterval of the range of integra- 
tion, or oscillates frequently across that range, 

This affliction can be diagnosed in many different ways. 

Sometimes a small change in an endpoint can render the 

integral improper, as in 

£ (lM1 in(u)du — 0.0989* ... -*£ ln|u]du = 1. 

Sometimes a small alteration of the integrand can render 
the integral improper, as in 

/_ dac/fx + I0~ lf1 ) = 314157.2654.^ -*/_ dx& = »« 

Sometimes the value of the integral is nearly independent 
of relatively huge variations in one or both of the end- 
points, as is J Q exp f-u 2 ] du «* VttV2 for all x > 10. Regard- 
less of the cause or diagnosis, nearly improper integrals 
are thfi bane of numerical integration programs, as we 
have seen. 

During the HP-34C 1 s design a suspicion arose that most 
integrals encountered in practice might be improper or 
nearly so. Precautions were taken. Now that experience 
has confirmed the suspicion, we are grateful for those pre- 
Cautions, They were; 

1. Avoid sampling the integrand at the ends of the range of 
integration. 

2, By precept and example in the Owner's Handbook, 
warn users against wild integrands, suggest how to recog- 
nize them, and illustrate how to tame them. 

The second precaution ignited controversy. Against it on 
one side stood fears that its warnings were excessive and 
might induce paranoia among potential customers. Who 
would buy calculator that he thinks gets wrong answers? 
Actually wrong answers were very rare, thanks in part to 
the first precaution, and many attempts to vindicate dire 
predictions about mischievous improper and nearly im- 
proper integrals were thwarted by unexpectedly correct 
answers like 

J*J In |nj du =a 0.999B ± 0,00021 

in l minutes al SCI 3, Or 

J exp(-u H } du = 0.886227 ± 0.0000008 

m 4 minutes at SCI 5, if the wages of sin be death, O Death, 
where is thy sting? 

On the other side stood a number of embarrassing ex- 
amples tike 
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i-IOO -, . 

J exp(-u"Jdu 

miscalculated as 0,0 ± 0.O0O0O0OOO5 in 14 seconds, 
Another, had we known it then, would have been Wood- 
yard 1 s example at the beginning of this article; I he correct 
answer 



0.12 -r 






u-1 



In u 



0,03649 ±0.0000000? 



in 23 minutes at FIX 7 differs from FIX 5 1 s wrung answer 
0,03062 in the worst way; the error is too small to be ob- 
vious and too large to ignore- Adding to the confusion 
were examples like 



AM 



"TV 



-2 In cos[ir] dun- = 1 - x 4 /60 + x a /480 + 



for which computation in SCI 4 produced ridiculous values 
like 4(0,1) = 0.95742 ± 0,00005, 4(0,01] = 0.58401 ± 
0.0000.!, and 4(0.001 1 = 0, all impossibly smaller than 1, 
This example ap pears to conde mn the /y key until the 
integrand f[u)=V -2 In cos[u 2 ]/u 2 is watched for small 
arguments u and seen to lose most of its figures to round- 
off, losing all of them for I u I ^0.003, despite an absence of 
subtract tuns that could be blamed for cancellation. Then 
the example appears to condemn the whole calculator. Who 
wants responsibility for a calculator that gets wrong 
answers? 

Don't panic [ The answers are wrong but the calculator 
is right, 

How to Tame a Wild Integral 

Forewarned is forearmed. Every experienced calculator 
user expects to encounter pathological examples like some 
of those above, and expects to cope with them. The ques- 
tion is not "whether" but "when"? And that is when atten- 
tion to detail by the calculator's designers is rewarded by 
the user's freedom from petty distractions t h Fit can only 
complicate a task already complicated enough. But like 
the dog that did not bark/ the absence of distracting de- 
tails may fail to be appreciated. That is why the examples 
explained below have been chosen — to illustrate the ad- 
vantages of liberated thought. Work them on your calcula- 
tor as you read them; don't skim them like a novel, Then 
you may come to think of your calculator the way I think 
of mine, as a trusted friend who stays with me when I 
need help. 

The integral A (x] above contains an integrand f(u) = 
V —2 in cos(u 2 l'u 2 that loses its figures when u becomes 
tiny. The problem Is caused by rounding cos(u ) to 1, 
which loses sight of how small u must have been. The 
solution compensates for roundoff by calculating f(u] 
as follows; 



Let y = uos u 2 rounded. 
If y = 1 then let f(u) = 1 
else let f(u) - V 



-2 In y.cos l y. 



The test for y = 1 adds four steps to the f- program and. 
provided In and cos -1 are implemented as accurately as 
on all recent HP calculators, the problem goes away. 

'SeeTrie las' *ew aaragrauhs nl The Sherlock Holmes story -Silver Blaze" by Conan Doyle 



V u/(u-1) - V\n{u) 



0.06 



ra> 




2w 2 /((w-1)(w+1)) - w/in(w) 



0.05 -r 



m 



0.05 




|r;: 



Fig. 10. Substituting w 2 foru turns the wsfd graph (a) into fne 
easy one (b). But do not replace {(w-1) (w+1)) by (w 2 -lj 
because roundoff errors introduce a spike, as shown in (c). 

Woodyard's example It has an integrand f(u) whose de- 
rative P{u) — k« as u — * 0. The graph of f(u) shown in Fig. 1 03 
looks like a lovers' leap. Stretching the u-axis near u=0 by 
substituting u = w z turns the precipice into the hummock 
shown in Fig. 10b and transforms the integral into an easy 
calculation: 



-f\ 2w! -S-) dw 

J ' iw-ij(w + l) lnw J 



The HP-34C computes this as 0.03649 ± 0,000005 in 100 
seconds at FIX 5 or 0.0364900 ± 0.00000008 in 200 seconds 
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at FIX 7. Do not replace (w-l)(w+l) by [w 2 -l) because 
the latter loses to roundoff half of its significant digits as 
w -* 1 and introduces a gratuitous spike into the integrand's 
Kraph shown in Fig. 10c. which was plotted on an HP-85. Do 
not worn' about w = or w = 1 because they don't happen, 
but do worry that as w — > 1 the integrand approaches tin 
reliable e ; -- - ~ j - == 0, This means that FIX 7 dis- 

plays about as many digits as could possibly be correct for 
all \v < 0.999, beyond which the JJ key draws few if any 
samples because it converges so fast. 

The graphs of expf — u"j over ^ u =£ 300 and of 
T/(u z + 10~ 10 )over-l ^u ^ 1 both resemble huddled mice 
with very long tails stretched out hundreds or thousands of 
times as Jong as their bodies. Plotting the graphs on a page 
of normal width is futile because the bodies get squashed 
into vertical whiskers. 

Most people who integrate such functions numerically 
cut off the tails. Thin tails can be cut almost indiscrimi- 
nately without much degrading the accuracy or the speed of 
integration, Such is the case for Jj exp f— u 2 )du, which J 
evaluates in less than, say. 4 minutes at SCI 5 provided that 
x, if bigger than 4 or 1 0, is cut back to someth ing between 4 
and 10, But J? x AvJ(u z + 10" 10 ) has too thick a tail to cut 
without losing accuracy or patience when x is large. That 
is why Procrustean methods are not recommended. Better 
to shrink the tail via an artful substitution like u = A + ^.tan v 
where a lies within the body of the mouse and fi Is roughly 
that body's width, Doing so with A = and ju = 1 changes 
J0 exp|-ir] du into 



-.■! i 

J-U 



expi tan^l |1 + tarry) dv 



which L evaluates in three minutes at SCI 5 even when 
x Is as big as 10 H \ Don't worry about tan tt/2 because it can't 
happen on a well-designed calculator. 

J\du/[ir + in "i 

benefits miraculously from the foregoing substitution 
when a = u and ft = 10 5 , but values near those do almost 
as well. 

Another technique might be called "subdivide and con- 
quer." It subdivides the range of integration into suhinter- 

1 1 upon each of which the integrand f(u) is tame, al- 
though i'(u) may l ook wild on the range as a whole. Fur 
example, f(u| a V' u 2 + 10~ 10 has a V-shaped graph prac- 
tically the same as that of lu . Kvahinling /„ 3 f(u) du 
accurately takes a long time if done with one press of 
f*< but subdividing the integral into 

/VhJidu +/jf[u)du 

takes two presses of/* and one of £+ but much less time, 
Subdivide and conquer works best when combined with 
apt substitutions. For example, if the formulas in Fig. 5 
were unavailable how would K(*) = f£ du/(l + u H4 l bq 
calculated? 

H'l J ( J du/(l + u b4 ) * |[" du/fl + u m ) ... subdivided 

= f t ] d«/|l + u 64 ] + J 1 w ff2 dw/(w w + l) ... u = l/w 



= f*{l * u fe2 )du|l - 

* 1 +/J (u &2 -u M J duf,i + u b4 l 

* 1 • -**) r'~' :: 

- 1.0O0401 7081 55 ± 1.2X10" 1 * 



merged via w-y 
. some algebra 

. u=v w to shrink 

a tail 



in 10 minutes at SCI 8. Thus we have calculated r 

4) to 13 significant decimals on a ten-sig- 
nificant-decimal calculator. 

Oscillatory integrals like Jp cosfln u] du sometimes suc- 
cumb to stretching substitutions like u ~ v 2 that damp the 
oscillations, but generally oscillatory integrals cannot be 
calculated accurately and quickly without sophisticated 
tricks beyond the scope of an article like this, A simple 
trick worth trying when the period of oscillation is known 
in advance is called folding, though it is really another 
instance of subdivide and conquer. Here is a didactic 
example. 



j 



W)0jt 



sin 2 u 



\ Li I VU + ff 



du = still running after over three 



hours at sci 5 



"INN 

t^O^ntr VU + Vu + JT 



if 

r- n ** ii 



$\n z v 



Vy-Hurr + S/v Mu- tt 



dv 



after being subdivided and with u « v + iitt. Exchanging 
J and ^ produces 



= f sin^ ■ 2! . 

Jfi n.n Vv+rnr + < 



Vv + lfff+ir 



dv. 



At this point a program should be written to calculate 
the sum. but because the example is didactic the sum 
collapses to yield 



, f * 6Q sinnj 



dv = 21,10204 ± 0,00007 



mow 



in 5 minutes a! SO 5 + 

Now for a final example drawn from life; 



-r 



du 



(a 2 + uj V[a J + u){b- + u][c 2 +ii) 



for a = 100,b=2,c = 1. 



This integral pertains to the electrostatic field about an 
ellipsoidal body with principal semiaxes a, b, c. 5 The 
ellipsoid is needle-shaped like an antenna or a probe, The 
classical approach transforms V into a standard form called 
an elliptic integral of the second kind and interpolates on 
two variables in published tables to gel a numerical value, 
The following approach takes less time. 

First transform the improper integral (Jq ) into a proper 
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one by substituting, say< u = (a 2 -c 2 ].'(l -v 2 ) - a 2 to get 

V = X / ' V[l-u-J/[u 2 + a) dv 
where 



k = 2'((a : '-(r) vV-b") = 2.00060U18 X 10 " ft 

/i = c'a - 0.01 

a = (tr-c^'ta^b 2 ) = |*0&t 200480 x in ' 
i 
Now, as always happens when a > > b > c, the integral 
is nearly improper because a and /i are both so nearly u. We 
suppress this near impropriety by finding an integral in 
closed form that sufficiently resembles the troublesome 
part of V r One candidate is 



W = k | dy/vV-H* = kin \v+Vv 2 +a) 



k In ((1+ V"l-htr),''[ J u + V'pt 2 + «J) 
B.401B1 880708 X 10 " fl 



Then 



V - W + X/ [V'(l ~v 2 \i[v 2 + a\ - i/Vt?*«) 4» 






= 7.78867525X10- S ± 1,3X10" l4 



after seven minutes at FIX 8. Don't worry about Vl-v 2 
as i/ — > 1 because the figures lost to roundoff are not needed 
and its infinite derivative doesn't bother the HP-34C. 

Conclusion 

A powerful mathematical idea has been placed at the 
disposal of people who will invoke it with fair confi- 
dence by pressing a button marked jy without having to 
understand any more about its internal workings than 
most motorists understand about automatic transmissions. 
Integrals thai might previously have challenged the 
numerical expert and a big computer now merely amuse 
the scientist or engineer, and tomorrow they will be rou- 
tine. And now those engineering students who do attend 
classes in numerical analysis need no longer be expected 
to memorize the names nor the remainder terms of quadra- 



ture formulas hut may instead be taught to use integra- 
tion wisely. 
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